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Abstract 



Weak boson fusion is a copious source of intermediate mass Higgs bosons at the LHC. 
The additional very energetic forward jets in these events are powerful background 
suppression tools. 1 analyze the decays if — >• 77 and H W^*^W^*^ — * e^ji^^T, with 
the latter a potential discovery channel, and the decay H t^t~ I'^h^'^T as a 
method for achieving the first direct measurement of a Higgs-fermion coupling. 

I perform parton level analyses of the signal and dominant backgrounds for each 
decay mode, and demonstrate kinematic cuts and other important tools neccessary to 
achieve an S/B > 1/1 rate in all cases. I also perform cross section calculations with 
additional gluon emission which provide an estimate of a minijet veto probability. 

I show that a 5a if —> 77 observation can be made for 110 GeV < Mh < 150 GeV 
with modest luminosity, order 40-50 fb~^ at low machine luminosity, overlapping the re- 
gion explored by the CERN LEP and Fermilab Tevatron. For 130 GeV < Mh < 200 GeV, 
I show that H W^*^W^*^ can achieve a 5a observation with S/B much greater than 
1/1 with extremely low luminosity, about 2-10 fb^^ over almost the entire range. This 
is the most promising search channel in the 130-200 GeV mass range. It overlaps the 
if — > 77 region and nicely complements the H — > measurement that can be 

made with very low luminosity in inclusive gg ^ H production. I further show that 
a Higgs-fermion coupling can be directly measured via the H ^ tt decay with only 
about 60 fb"^ (assuming low luminosity running). 
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Chapter 1 
Introduction 

1.1 Where Are We? 

The year 1999 finds the particle physics community on the cusp of understanding the 
fundamental particles and forces of nature. We have, over the past three decades, 
formulated a theory governing the electromagnetic, weak and strong forces and their 
constituents, called the Standard Model, which now predicts with uncanny accuracy 
the results of most particle physics experiments, and there is a considerable volume of 
data with which to compare [1]. During this time numerous observations have been 
made that suggested new physics beyond the Standard Model, or that it was wrong, 
but all have essentially vanished as statistical fiuctuations or experimental fiaws. 

Yet while the community continually rejoices in the predictive success of the Stan- 
dard Model, it also faces the somber realization that the theory is unsatisfactory, that 
it is incomplete at a fundamental level. Aside from the deeper questions of the large 
number of input parameters required to define the theory (of order 20) [2], the intima- 
tion of unification of forces at energy scales higher than can be probed [3] , and perhaps 
even the necessity of such principles as renormalizability, there are the more immediate 
problems of electroweak symmetry breaking (ESB) and mass generation, on which this 
dissertation focuses. 
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The Standard Model hinges on an explanation of these observations at a time when 
the community is on the verge of having the technical capability of exploring some of 
these issues with conclusion. It remains to be seen whether or not the Higgs mechanism, 
the (currently) most viable explanation of ESB, is correct. It is certain that this area has 
not yet been fully explored phenomenologically, and this dissertation outlines strategies 
for experimental searches designed to lay the issue to rest with efficiency and certainty. 

In this Chapter, I introduce the physics of the electroweak sector of the Standard 
Model which describes all presently known interactions and forces other than gravity, 
explaining its structure and the mathematical limitations in formulating the theory with 
only the particles observed in nature. I describe how the observed ESB of the theory 
can be obtained via addition of a simple Higgs sector, and show the Standard Model 
couplings that result. I explain the origin of mass in the theory, for both known particles 
and the theoretically introduced but unobserved Higgs scalar, and why a light Higgs 
mass is expected from radiative corrections to the masses of the weak gauge bosons 
in light of the current status of electroweak (EW) precision fits to data. As a prelude 
to introduction of the minijet veto tool in Chapter 2, I discuss the basics of QCD, the 
strong force responsible for the interactions of colored particles. Finally, to introduce my 
calculations for a Higgs search in weak boson fusion I give a brief description of parton- 
level Monte Carlo calculation of cross sections, and an introduction to madgraph [4], 
a tool I use for generating matrix elements in FORTRAN. 

1.2 The Standard Model 

The Standard Model (SM) has evolved from a field theoretical description of only elec- 
tromagnetism (Quantum Electro-Dynamics, or QED) [5], to a much richer theory en- 
compassing also the weak force, first observed in nuclear beta decay, and the strong 
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force, which binds together nucleons and their constituents, quarks, which are dis- 
cussed further in Section 1.5. All known particles in nature (other than the graviton) 
are described by the SM. Two classes of fermions, leptons and quarks, make up matter. 
Bosonic particles carry the three forces which act between them: the massless photon, 
7, for electromagnetism; the neutral Z and charged VF^, both of which are massive, for 
the weak force; and the massless gluons, gi, for the strong force. 

I begin by describing the group theoretic structure of the SM - the symmetries 
SU (3)c ® SU (2)i (8) [/ (l)y , where c refers to the color group of QCD, L refers to a left- 
handed group structure, and Y is the hypercharge group [6,7]. To this gauge structure 
one adds a scalar SU{2) doublet of hypercharge -|-1, called the Higgs doublet, which 
breaks SU{2)l <8) U{1)y — U{1)em [7,8]: massive weak gauge bosons W^,Z plus 
the electromagnetism of massless photons. The left-handed nature of the SU (2) group 
comes from direct observation of left-handed doublets of matter fields in nature, the 
quarks and leptons, along with right-handed singlets: 




dR ur 



gr 



V e / 




Sr Cr 



fJ'R 




bR tR 



Tr 



[l.l) 



While all matter fields transform non-trivially under [/(l)y, only the left-handed dou- 
blets transform non-trivially under SU{2)l, and only quarks and gluons carry the con- 
served charges of QCD. Note that while we have observed all ^ the particles of the three 
complete generations shown above, only the first generation is found in normal matter: 
particles of the second and third generations are more massive and unstable, decaying 
in all cases to particles of the first generation. The SM fails to predict the masses of all 
matter fields. 



^The tau neutrino, has as yet not been observed directly, but is believed to exist based on 
indirect evidence, such as the total decay rate r{Z — > uu). 
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Each symmetry group has its own generators of transformation, and thus its own 
force-carrying bosonic particles, in the adjoint representation. For the color symmetry 
of QCD there is an octet of bosons, for SU{2)l there is a triplet, and for U{1)y there 
can be only one. One may write a Lagrangian for the SM gauge and matter fields as [9] 

C'SM — + + ^GF + ^FP , (1-2) 

where jCq is the kinetic term for gauge fields, £/ is the kinetic term for fermions, which 
also (via minimal coupling) describes the interactions of fermions with gauge bosons, 
and Cgf and Cpp ai'c gauge-fixing pieces I do not address here. A fundamental principle 
in the construction of this Lagrangian is that of gauge invariance - the only theories of 
vector bosons that can be renormalized are gauge theories. Thus, the gauge principle 
is a basis for the predictive power of the SM. 

The first problem with this Lagrangian is that by construction all gauge and matter 
fields must be massless. Inserting mass terms by hand leads to disaster, as gauge boson 
mass terms are not gauge invariant, and matter field mass terms will not respect the 
SU{2)l ® U{1)y symmetry. However, in nature we observe all matter fields to have a 
mass, as well as the weak bosons W^, Z. 

The second problem is that of bad behavior at high energy for massive W and Z 
bosons. This same feature was present in the Fermi theory of four-fermion interactions, 
the predecessor to the SM [10]. Weak processes, first observed in nuclear beta decay, 
were described via the coupling Gp and an expansion in powers of the center-of-mass 
(CoM) energy, E^, for the scattering amplitude. Gp was measured in experiment, and 
the theory worked very well, but was doomed to ultimate failure as it predicted violation 
of unitarity for E ^ 600 GeV. To regulate this behavior, a modification to the theory 
was proposed, whereby the interaction of four fermions is replaced by a force-carrying 
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massive boson, a charged W^, which predicts an amphtude of the form 



where g is the couphng, related to the Fermi couphng by = J^^m^- inspection we 
see that the low-energy {E <^ M]v) behavior asymptotically approaches that described 
by the Fermi theory, amplitude proportional to E"^, while at high energy, the amplitude 
approaches a constant, thus avoiding unitarity violation. High-energy data supported 
this description to incredible accuracy, and eventually the massive bosons were 
directly observed along with a massive neutral Z boson, which confirmed an SU (2) l ® 
U{1)y structure as integral to the SM. 

This theoretical fix to Fermi theory merely put off the problem, however: calculation 
of SM weak boson scattering at high energies, e.g. W'^W~ W'^W~, reveals that 
unitarity is again violated at CoM energies of order ~ 1 TeV. The theoretical success 
of the Higgs mechanism lies in supplying not only masses for gauge bosons and matter 
fields, but also in regulating this bad behavior in weak boson scattering. Including 
Higgs exchange in W~^W~ — > W~^W~ scattering is exactly what restores unitarity. 



1.3 Elect roweak Symmetry Breaking 

To obtain massive gauge bosons the SU{2)l ®U{1)y group structure must be broken. 
We must obtain after the breaking three massive gauge bosons with the observed mass 
ratio, and a massless photon, reflecting a leftover, unbroken U{1)em- Massless vector 
fields have two degrees of freedom, whereas massive fields have three, so one must 
introduce additional degrees of freedom that account for the third helicity state of 
the bosons that become massive. Following mostly the notation of Ref. [7], this can 
be done most simply by introducing an additional spin-0 SU{2)l complex doublet of 
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hypercharge Y — +1, written as 

The electric charges of the components are known via the relation Q — T3 + ^, where T3 
is the third component of the SU{2)l quantum number, ±|, and Y is the hypercharge 
of the doublet. This doublet "spontaneously" breaks the symmetry when it acquires 
a non-zero vacuum expectation value. It is inserted into the theory via the additional 
Lagrangian term 

Ch = (1.5) 

V{^) = //2($t$) + A($t$)2 , (1.6) 

with minimal coupling to the gauge fields via the covariant derivative 

D,^d, + ig^W; + ig'^B, {i = 1, 2, 3) , (1.7) 

where FF^ are the massless SU{2)l gauge fields and is the massless hypercharge 
gauge field. This potential is the simplest that can be written which gives a non-zero 
expectation value for the field and is bounded from below, provided that //^ < and 
A > 0. It has a global minimum at 

\{wl + wl + <t>' + wl) = (1.8) 

One may choose which component of $ lives at this minimum, and for convenience 
choose it to be 0: 

Wi = W2 = W3 = , 0^ = — - = . (1.9) 

A 

This is the vacuum expectation value (vev) of the Higgs field. $ can now be written as 
small perturbations about this minimum. 



,.exp|!^^:)4(^^°^|. (1.10) 
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Since $ and the component fields respect the SU{2)l symmetry, one may make a 
transformation to the Unitary gauge and write the simpler expression 

Insertion of this resulting expression for $ into the kinetic term in Lh, we find the 
Lagrangian terms 



1 



1 



Ch = ... + -{v + Hfg' [W;^!^ + W^^W^f^l +-{v + H)\g'B^ - gW^^){g'B^^ - gW'^ ■ 

(1.12) 

Defining the observed fields as — -^{W^ ^ W^), we immediately find a mass term 

{■^v'^g'^W^W''^), a trilinear coupling term (^vg'^HWj^W^^^), and a quartic coupling 
term {^g'^HHW^W~'^). The W boson mass may be immediately read off to be Mw — 
\vg. 

The part of the second term's expansion may be expressed as a mass mixing 
matrix 

/ ^ 1^ 1^3/. ^1 
-gg' g'^ )\b^ 

which has eigenvalues and \\/ g^ + g''^ , the observed masses of the photon and Z 
boson. Their fields are similarly defined in terms of and via a mixing angle, 
called the Weinberg angle, given by the relation 



1 2,T./3 ^^19 -99 



^Awi,B,) 



(1.13) 



V9 + 9 



This leads to the relation for electric charge, 

e — g sin 9w — 9' cos 9w ■ (l-l^) 

There are similar terms for HZ^Z^ and HHZ^Z^ couplings. Note that there are no 
trilinear HA^j^Ai^ coupling or HHA^j^Ai^ quartic couplings! Thus, the i7 — > 77 decay is 
a rare process induced by loop diagrams. 
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The Higgs mechanism also predicts a relation between the W and Z masses, which 
is a crucial test of the SM: 



= 1. (1.16) 



~ Ml cos2 ew 

Another crucial test of the SM is the requirement to observe the trilinear HVV couplings 
- only models with scalars that acquire a vev can produce this feature. Thus, observation 
of a trihnear (j)VV coupling would identify the Higgs as the agent responsible for ESB. 

It is worth noting that the four additional fields added via $ have been rearranged 
such that three of the fields have disappeared while three degrees of freedom have shown 
up as the longitudinal components of the newly-massive weak gauge bosons, and the 
fourth field has become a real, interacting particle. To see how the Higgs field itself 
acquires a mass, Eq. (1.11) must be inserted into the potential V{^) of Ch to obtain 
additional terms in the Lagrangian: 

^X{vH + H^)^ (1.17) 

which give Higgs self-interactions {HHH and HHHH couplings) as well as a Higgs 
mass, niH — v\/2\. It is very important to realize that the only places in the theory 
where A appears are in the Higgs self-couplings and the Higgs mass term. Thus, ex- 
periment cannot determine the Higgs mass directly without observing the Higgs or a 
process involving the Higgs propagator. Other methods can, however, give us clues as 
to its mass, as explained in the next Section. 

Fermion masses, while distinct from the issue of ESB, can also be generated from a 
Higgs SU{2) doublet via Lagrangian terms of the form Y^Fl^cIr ^^(^ ^fFL^fk-i ^-^d 
their Hermitian conjugates, where are left-handed fermion doublets, are right- 
handed fermion singlets, and Yf are the Yukawa couplings for up-type and down-type 
fermions, respectively. Taking the first-generation leptons as an example, one may write 
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the Lagrangian terms 

I 



/ ^ 



V2 



eR -YeeR^{0,v + H) 



e 



:i.i8) 



\v + H ) 

= ... -^{v + H)eLeR -^{v + H)eReL, 

which yield a mass for the electron, ■^Y^.v, and an electron-Higgs coupling, naturally 
proportional to the mass, Each fermion has a different Yukawa coupling Yf, and 
for the quarks the values are different for Y^ and Y^. The Standard Model does not 
explain why the Yukawa couplings exhibit the seemingly unrelated values that they do; 
here, the mechanism merely provides mathematical consistency for construction of the 
theory. It is quite unsatisfactory in explaining why the fermions have their observed 
masses. 

It is also possible to introduce two Higgs doublets to the theory. Typically, this is 
done such that each Higgs acquires a vev. By construction this must give the same ESB 
and gauge boson masses, but one doublet's vev will give masses to up-type fermions 
and the other's to down-type fermions. While there has been much research on general 
two-Higgs doublet models, this scenario is more commonly explored in the context of 
the minimal supersymmetric Standard Model (MSSM). I discuss this model briefly in 
Appendix F, and also extend the SM results of Chapter 5 to the MSSM case. 



1.4 The Higgs Mass 

While the Standard Model does not predict the Higgs mass, there are theoretical con- 
straints on the allowed mass range. First and foremost is the unitarity limit, obtained 
by examining weak boson scattering at high energies. Unitarity is maintained if the 
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Higgs mass is less than about a TeV [11], 

ruH^l TeV (unitarity) . (1.19) 

A closer examination of the HHHH coupling suggests a stricter limit [12,13]. This 

coupling runs according to the relation 

1 1 3 , /A\ 

+ — log - , (1.20) 



A(At) A(A) 27r2 \pL 
where /i is some low scale at which we measure A(//), and A is a higher scale. The first 
requirement is that A(A) never be negative, an unphysical value. This leads, after some 
algebra, to an upper bound on the Higgs mass. 



Mh < 



^31og(-^V 



;i.2i) 



Mh) 

This is known as the "triviahty" condition, because if the SM Higgs sector is to be vahd 
for arbitrarily large energies (efi:ectively A ^ oo with A (A) finite) then \{Mh) — 0, 
making it a trivial theory. For example, for A = 10^^ GeV, one finds the limit Mh ^ 
175 GeV. For small values of A there is still a limit, although in this case the running 
of A(/i) enters a nonperturbative region and the above equation becomes significantly 
modified; lattice simulations of the small-A limit estimate that mn ^ 630 GeV = A [13], 
much less restrictive than for A — > oo. 

The requirement of vacuum stability places instead a lower bound on the Higgs 
mass [12]. Here, one-loop corrections to the Higgs potential, which must have an abso- 

i ° 

lute minimum at = , yield the relation 

M'h > ^^^(2M^ + M| - mt) log (^) > -4500 log (^) GeV^ . (1.22) 

That this bound is negative and therefore not a real constraint is a result of the fact 
that rrit ^ Mz- Two-loop corrections, however, are sizeable and yield a positive-definite 
lower bound. 
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An important feature to note of the above relations is that for small A, there is 
essentially no bound on Mh- Thus, while there are theoretical arguments for expecting 
A to be large, e.g. the GUT or Planck scale, one cannot take these expressions as real 
restrictions since there is no experimental confirmation of a desert. 

Instead, the most important hints of what the Higgs mass may be come from in- 
direct measurements in present-day experiments - precision fits to electroweak observ- 
ables. Both the Higgs and the top quark appear in loop radiative corrections to the 
gauge boson propagators. Now that the top quark has been observed and its mass 
measured, precision fits can place upper and lower confidence limits on the Higgs mass. 
Unfortunately, while these radiative corrections are a function of m^, and therefore very 
sensitive to the top quark mass, they are a function of only the logarithm of the Higgs 
mass, log(mi^), making the sensitivity very weak. Nevertheless, these fits give us an 
important clue: the Higgs mass should be low, on the order of 100 GeV, as opposed to 
several hundred GeV or close to the unitarity limit. Recent (1999) fits give the best-fit 
value for the Higgs mass to be [12] 

mjj = 1071^^ GeV , (1.23) 

with a 95% C.L. upper limit of niH ^ 255 GeV. While the uncertainty on this is quite 
large, the point is that indirect data suggests the Higgs has a fairly low mass, at the 
low end of the "intermediate" range (100-200 GeV). This should not be surprising, as 
from the historical example of Fermi theory one would expect the regulating physics 
to step in at a scale much lower than the unitarity limit; the weak bosons have masses 
mw = 80.35 GcV and mz = 91.19 GeV, much lower than the upper bound of 600 GeV 
imposed by unitarity. 

Direct searches have established an ever-increasing lower limit on the Higgs mass 
as machine energies have risen. The current bound, established by the LEP collider at 
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CERN in 1998, places the 95% C.L. limit [14] 



ruH > 99 GeV . 



(1.24) 



In its next (and last) two runs through 2000, LEP will either discover a Higgs up to a 
mass of about 105 GeV, or place a 95% C.L. exclusion limit of about 110 GeV. 



Quantum Chromo-Dynamics, or QCD, accounts for the remaining sector of the SM. 
It is an unbroken local gauge symmetry, described by the non-Abelian SU{3)c gauge 
group, where c refers to color [15]. The non-Abelian aspect means the generators 
of SU{3) transformation do not commute. This leads to some interesting properties 
which will be relevant for an investigation of the minijet veto tool and calculations I 
perform for QCD processes. The non-Abehan nature leads to gauge bosons that carry 
the conserved charges. Thus, the eight massless gluons of SU{3)c QCD are also color- 
charged particles. In contrast, the massless gauge boson of QED, the photon, carries 
no electric charge. Since interactions are governed by conserved charges, photons may 
not have self-interactions, but gluons must. 

The value of the conserved charge in a gauge theory is modified by higher-order 
corrections to the vertex, as well as vacuum polarizations in the gauge boson propagator. 
In QED, these can only be fermion loops, which lead to a running of the coupling 
constant. 



where is the energy scale being probed and is some reference scale where q;(/x^) is 
finite and known. It is easy to see that at low energy, the QED coupling is essentially 
constant, but at very high energy it becomes large - calculations in this regime become 
nonperturbative. It turns out that this regime is well beyond the reach of experiment. 



1.5 QCD: Theory and Calculation 




(1.25) 
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In QCD, the opposite happens. Because gluons self-couple, gluon loops in the vac- 
uum polarization propagator corrections must also be calculated. This leads to a run- 
ning coupling of the form 



Here, the coupling blows up as becomes small, but approaches zero as — > oo; 
this feature, which is general of non-Abelian gauge theories, is known as "asymptotic 
freedom". Because of asymptotic freedom, QCD calculations instead become nonper- 
turbative at low energy, which are referred to as the soft regime or soft effects. Our 
minijet veto tool naturally requires soft jets, which may be nonperturbative in the ex- 
perimental phase space regions of interest. I explain how this can be dealt with in the 
next Chapter and in Appendix C. 

The self-couphng of gluons has an additional effect on colored particles: lines of flux 
between color charges attract each other and so become constricted, such that the lines 
are confined to a narrow "fiux tube" connecting the two physical particles. In QCD, the 
energy density per unit length in a flux tube is approximately constant, leading to a color 
fleld potential that grows proportional to the separation between the colored particles, 
V{r) ~ r, and ultimately results in confinement. That is, colored particles must exist 
only in "colorless" bound states, such as color/anti-color mesons or red/ green/ blue (or 
anti-red/anti-green/anti-blue) baryons, such as the proton and neutron. High energy 
collisions may separate colored particles; if the momentum transfer is sufficient, we may 
calculate this hard scattering process perturbatively. Additional quark/gluon splittings 
may then occur, called parton showering because the relatively large probability of 
this occurring typically results in many more final state particles than the 2 — > n hard 
process under consideration. Showering, too, may be calculated perturbatively, down to 
momentum transfer values of order ~ 1 GeV^. Below this scale, the colored partons 
separate from each other sufficiently that one must consider the long-range behavior of 
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QCD, which is inherently non-perturbative. The kinetic energy of the separating quarks 
and gluons is converted into potential energy in the flux tubes maintaining the color 
connections. At large enough distances, when the energy in the flux tube is suflicient, a 
pair of colored particles is created out of the vacuum, which then separate themselves. 
This continues until not enough energy is left in a flux tube to create new quark/anti- 
quark pairs and the system condenses to a collection of colorless bound states. This 
process is known as hadronization. 

Showering and hadronization lead to important consequences for collider physics. 
For collisions involving either initial- or final-state colored particles, for example uu — > 
d'd', these processes lead not to the final state with only two observed particle tracks, 
but to two clusters of hadrons, which may be grouped together into jets. Showering 
and hadronization are well-understood phenomenologically, but for practical purposes 
can be simulated only via tuned Monte Carlo programs [16] (Monte Carlo refers to 
the random nature of the numbers used in generating phase space points). Various 
programs such as phytia [17], HERWIG [18], and ISAJET [19], exist that take into 
account successively the hard scattering, parton showering, and hadronization. They 
can also simulate showering for collider processes involving more than two final-state 
particles in the hard scattering process. For the background processes I examine in 
the following Chapters, however, these programs would first require additional tuning 
because the acceptance requirements imposed on the final states do not necessarily 
satisfy the coUinearity assumptions of the programs' approximations. For example, 
PYTHIA and HERWIG are known to give much too small a cross section for Z+ hard 
jets [20]. Thus, input from parton-level studies with full tree-level matrix elements is 
required, which this dissertation provides. 

While the calculations are more reliable when using tree level matrix elements for 
a hard scattering process, one is then forced to make an approximation in place of full 
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parton shower simulation. For the simulations here I equate final-state colored particles, 
which arc required to be well-separated, with final-state jets of a given radius in the 
lego plot ^ centered around the single particle. This parton-level-only approximation is 
sufficient for proof-of-concept of search strategies and for developing tools to enhance 
a signal relative to its backgrounds. More sophisticated parton showering calculations 
with complete matrix elements for the hard scattering processes will eventually be 
required for comparison of theory with experimental data. 



1.6 Monte Carlo Techniques 

Calculation of a cross section in field theory is an integration of a squared matrix element 
for a given process, summed over final-state quantum numbers such as helicity or color, 
averaged over the same initial-state quantum numbers if the initial-state is unpolarized, 
and integrated over the phase space of the initial- and final-state particles [21]. A total 
cross section is the sum of cross sections for all distinct initial/final-state subprocceses. 
This may be written as 

(7 = J dxidx2 Yl fai{xi)fa^{x2)^ J d^nQicuts)Y^\M\'^{subproc) , (1.27) 

subproc 

where the summed/averaged/squared matrix element is 

4 #COlors(aia2j colors pol'zn 

(the I comes from both initial-state particles ai and 02 being partons with two polar- 
ization degrees of freedom), the Lorentz-invariant phase space (LIPS) is 

d^^{P-p,...p^) = f[ [j^^^ (27r)V(P - , (1-29) 



^Rectilinear plot of the convenient detector variables pseudorapidity (r/) v. azimuthal angle {(f)). 
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the faii^i) are the parton distribution functions, which give the probabihty to find 
parton with momentum fraction Xi inside the incoming hadron, and s = s ■ Xi ■ X2- 

For hadron coUiders, such as the LHC, it is impossible to calculate any cross section 
analytically, even for total rates, because the squared amplitude must be folded with 
structure functions for the incoming (anti-) protons. In addition, the calculational goal is 
typically not a total rate, but instead an integration over only the part of the phase space 
where the final-state particles are visible in the detector. Further cuts on phase space 
may also be enforced to extract a signal from the total background. As such, hadron 
collider calculations require the use of a computer to perform the integration over a very 
complicated phase space configuration. We may approximate the integration over the 
LIPS by transforming to a more convenient set of variables, preferably across which the 
differential cross section is fairly flat, and sampling this phase space with quasi-random 
numbers. This transformation involves a Jacobian that must be reevaluated at every 
phase space point: 

3n-2 

—dxidx2d^n ^ J Yl dvi. (1.30) 
2* i=i 

We then instead evaluate the expression 

^-^E-^QMs) E fa,ix,)fa,{x2)Y,\M\\subproc), (1.31) 

{n} subproc 

where the matrix element squared is evaluated numerically at each phase space point. 
Typically, the phase space variables dpi^dpiydpi^... are very inefficient, and an expres- 
sion for the Jacobian transform to detector variables such as pt, t], is performed by 
hand before coding, pt is simply a convenient variable for detector capability, and r] 
is useful because it expands the angular region where small angle scattering results in 
a divergent cross section, thus making it more efficient for the integration routine to 
avoid. The goal is to choose a set of phase space variables such that the cross sec- 
tion is nearly flat across < n < 1 for each variable. In practice, this is difficult to 
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achieve, but a transformation to detector variables is vastly more efficient than using 
rectilinear variables. 1 outline the general structure of our Monte Carlo programs and 
our technique for generating matrix elements in Appendix E. 
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Chapter 2 
The Weak Boson Fusion Signature 

2.1 Introduction 

For the intermediate mass range, Mh f« 100 — 200 GeV, most of the hterature has 
focusscd on Higgs production via gluon fusion [23] and tiH [24] or WH[ZH) [25] asso- 
ciated production. Cross sections for Higgs production at the LHC are well-known [23], 
and while production via gluon fusion has the largest cross section by almost one order 
of magnitude, there are substantial QCD backgrounds but few handles to distinguish 
them from the signal. Essentially, only the decay products' transverse momentum and 
the resonance in their invariant mass distribution (if it can be constructed) can be used. 
The second largest production cross section for the standard model (SM) Higgs boson 
is predicted for weak-boson fusion (WBF), qq qqVV QqH, shown in Fig. 2.1. 
WBF events contain additional information in the observable quark jets. Techniques 
like forward jet tagging [26-28] can then be exploited to reduce the backgrounds. I 
discuss the basic signal process and characteristics, including the results of the forward 
jet tagging technique, in Sec. 2.2. 

It is necessary to study other production channels for several reasons. For instance, 
electroweak symmetry breaking (ESB) and fermion mass generation may be less inti- 
mately connected than in the Standard Model (SM) and the coupling of the lightest 
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Figure 2.1: Weak boson fusion Higgs production. 

Higgs resonance to fermions might be severely suppressed. In this case, neither gg ^ H 
fusion nor ttH associated production would be observed. Once the Higgs is observed 
in both gg ^ H and the WBF process qq qqH, where the Higgs is radiated off vir- 
tual Ws or Z's, the cross section ratio of these modes measures the ratio of the Higgs 
couphng to the top quark and to W,Z. This value is fixed in the SM, but deviations are 
expected in more general models, like super symmetry with its two Higgs doublets [29]. 
As I shall demonstrate, the WBF channel may yield a quicker discovery, requiring only 
2-30 fb~^, depending on the Higgs mass, which is less than the integrated luminosity 
needed for discovery in the gg ^ H channel in some decay modes. Since the Higgs 
can be discovered in the WBF channel at relatively low integrated luminosity, a fairly 
precise measurement of various Higgs couplings can be obtained with 100 fb~^ of data 
or more. 

Another feature of the WBF signal is the lack of color exchange between the initial- 
state quarks. Color coherence between initial- and final-state gluon bremsstrahlung 
leads to suppressed hadron production in the central region, between the two tagging- 
jet candidates of the signal [31]. This is in contrast to most background processes, which 
typically involve color exchange in the t-channel and thus lead to enhanced hadronic 
activity in the central region. We may exploit these features, via a veto on additional 
soft jet activity in the central region [32]. I describe the idea of a minijet veto in Sec. 2.3, 
including different techniques for estimating veto probabilities. 



20 



2.2 The WBF Signature 

The signal can be described, at lowest order, by two single-Feynman-diagram processes, 
qq qqiWW, ZZ) — >■ qqH , i.e. WW and ZZ fusion where the weak bosons arc emitted 
from the incoming quarks [33], as shown in Fig. 2.1. As the Higgs is a spin-0 particle, 
its decay may be treated separately from production, either as a branching ratio or via 
a decay matrix element which is squared and summed over helicities separately, then 
simply multiplied to the production cross section. Thus, we can discuss general features 
of WBF Higgs production independent of the decay channel. ^ 

The first task is to identify the search region for these hard Hjj events. Features of 
the signal are a centrally produced Higgs which tends to yield central decay products, 
and two jets which enter the detector at large rapidity compared to the decay products. 
I start out by discussing three levels of general cuts on the qq — > qqH signal, before 
considering decay products and their identification. This procedure makes exphcit the 
source of the major signal reduction factors encountered. 

The basic acceptance requirements must ensure that the two jets and the Higgs 
decay products are observed inside the detector. The outgoing quarks/gluons must lie 
within the reach of the hadronic calorimeter, with sufficient pt to be identified as jets. 
As the decay products are typically leptons or photons, they are required to lie within 
the electromagnetic calorimeter. In addition, the jets and observable decay products 
(generically denoted by X) must be well-separated from each other: 

Pt, > 20 GeV , \r]j\ < 5.0 , A% > 0.7 , 

|77x|<2.5, ARjx>0.7. (2.1) 

^Physical parameters, parton distribution functions and the factorization scale are chosen as in 

Appendix A. Detector simulation is included to the extent of energy /momentum resolution of the 
final-state particles via Gaussian smearing prescriptions based on the LHC detector expectations; see 
Appendix B for details. 
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Figure 2.2: Normalized minimum and maximum distributions of the two jets in 
signal qq qqH, i7 — 77 events at the LHC, with Mh = 120 GeV. No cuts are 
imposed; the vertical dashed lines represent the minimum detector acceptance cut 
of Eq. 2.1. 

The maximum and minimum jet pt and r] distributions of the Hjj signal in Figs. 2.2 
and 2.3 demonstrate that these basic acceptance requirements remove little of the signal. 
The figures are shown for Mh — 120 GeV, but this is true regardless of mass, which 
is not surprising as the pt's of the jets are governed by the scale of the weak boson 
masses and are thus rarely small. Similarly, the maximum rapidity distribution of the 
Higgs decay products, shown for the case of photons for simphcity in Fig. 2.4, tend to 
be very central, well within the reach of the electromagnetic calorimeters, which will 
extend to about 2.5 units in 77 for both CMS and ATLAS. The QCD, backgrounds, 
on the other hand, seen in the same figures, typically have more central jets and more 
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minjets I'? I max.jets 

Figure 2.3: Normalized pseudo-rapidity distributions of (a) the most central and (b) 
the most forward of the two tagging jets in •yyjj events at the LHC. The generic 
acceptance cuts of Eq. (2.1) are imposed. Results are shown for the qq — > qqH signal 
at niH — 120 GeV (solid hne) and for the irreducible QCD 77 background (dashed 
hne) . 

forward photons. Slightly more than half of all signal events pass these basic cuts, as 
seen in the first two columns of Table 2.1. 

The fact that irreducible QCD backgrounds typically give higher-rapidity X's as 
compared to the Higgs signal is due to the bremsstrahlung nature oi W, Z and 7 
emission, which occur at small angles with respect to the parent quarks, producing 
7's or leptons forward of the jets. To illustrate this more clearly, wc may define a 
new variable, a shifted photon rapidity, r}*, which is the rapidity of the photons with 
respect to the center of the tagging jets, 77*. = 77-y. — ^{rjtagi + Vtag2), which is shown for 
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Figure 2.4: Normalized pseudo-rapidity distributions of the most forward of the two 
photons in 77 events at the LHC. The generic acceptance cuts of Eq. (2.1) are im- 
posed. Results are shown for the qq — > qqH signal at tuh = 120 GeV (solid line) and 
for the irreducible QCD 77 background (dashed line) . 

Hjj, H —>■ jj signal and QCD ■yyjj background events with the basic cuts of Eq. (2.1) 
in Fig. 2.5. An obvious cut is to reject events that typically have Higgs decay products 
outside the central region defined by the tagging jets. Thus, at the second level of cuts 
both X's are required to lie between the jets with a minimum separation from each jet in 
pseudorapidity Arjjx > 0.7; and the jets are required to occupy opposite hemispheres: 

+ 0.7 < 77x1,2 < 'nj,max - 0.7 , Tjj^ ■ Tjj^ < . (2.2) 

Finally, Fig. 2.6, which shows the tagging jet separation in 77 for Hjj, H ^ 'j'j and 
QCD 77JJ events, reveals that the separation is drastically wider for the signal. The 
narrow separation of the tagging jets is a general feature of QCD processes at the LHC, 
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Figure 2.5: Normalized pseudo-rapidity distributions of both photon rapidities with 
respect to the center of the tagging jets (as discussed in the text) in events at 

the LHC. The acceptance cuts of Eq. (2.1) arc imposed. Results are shown for the 
qq qqH signal at mn = 120 GeV (solid line) and for the irreducible QCD 77jj 
background (dashed line). 

which are largely gluon-initiated; large separations require higher invariant masses, and 
thus restrict the processes to the region of Feynman x where quarks dominate. Thus, 
at the third level of cuts, which is also the starting point for our consideration of the 
various backgrounds, a wide separation in pseudorapidity is required between the two 
forward tagging jets, 

A77ta3, = |r/j, -r^,,! >4.4, (2.3) 
leaving a gap of at least 3 units of pseudorapidity in which the X's can be observed. 
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Figure 2.6: Normalized pseudo-rapidity distributions of the width of the allowable Higgs 
decay product observable gap region 77, Irjtagi — Vtag2\ — 1.4, in 77 jj events at the LHC. 
The acceptance cuts of Eqs. (2.1,2.2) are imposed. Results are shown for the qq qqH 
signal at ttih = 120 GeV (solid line) and for the irreducible QCD 77JJ background 
(dashed hue). 

This technique to separate weak boson scattering from various backgrounds is well- 
estabhshed [26-28,32,34-36], in particular for heavy Higgs boson searches. Table 2.1 
shows the effect of these cuts on the signal for a SM Higgs boson over the mass range 
niH = 100 — 200 GeV. Overall, about 29% of all Hjj events generated in weak boson 
fusion are accepted by the cuts of Eqs. (2.1-2.3). These cuts form the core identification 
technique for a WBF Higgs signal. Additional, unique cuts used to distinguish the very 
different cases X = 7, PV, r will be discussed in the appropriate Chapters. 
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Table 2.1: Signal inclusive cross sections (pb) for Hjj events of various Higgs masses 
in pp collisions at ^/s — 14 TeV. Results are given for no cuts and successive cuts of 
Eqs. (2.1-2.3). For the angular cuts of Eq. (2.1) on the Higgs decay products, H is 

assumed to decay into two massless particles {bb, t^t^ , 77, etc.). 



Higgs mass 
(GeV) 


inclusive 


+ Eq. (2.1) 


+ Eq. (2.2) 


+ Eq. (2.3) 


100 


4.8 


2.7 


1.62 


1.21 


120 


4.1 


2.4 


1.45 


1.10 


140 


3.5 


2.1 


1.30 


1.00 


160 


3.1 


1.8 


1.18 


0.92 


180 


2.7 


1.6 


1.07 


0.84 


200 


2.4 


1.4 


0.97 


0.78 



I should note that this calculation of the Hjj signal process has very low uncertain- 
ties, as reflected by its stability against factorization scale uncertainty: variation of the 
factorization scale by a factor of two changes the 2- jet cross section at the level of for- 
ward tagging cuts (Eqs. (2.1-2.3)), shown in the last column of Table 2.1, by < ±10%. 
This implies that next-to-leading order (NLO) corrections to this process are expected 
to be small, which is borne out by calculations found in Ref . [30] . 

2.3 Distribution Patterns of Additional Radiation 

A further characteristic of EW vs. QCD scattering can be exploited, namely the absence 
of color exchange between the two scattering quarks in the qq — > qqH signal process, 
^-channel color singlet exchange in the EW case leads to soft gluon emission mainly 



27 



in the very forward and very backward directions, whereas QCD processes are domi- 
nated by t-channcl color octet exchange which resuhs in soft gluon radiation mainly in 
the central detector. It was hoped that resulting rapidity gaps in signal events (large 
regions in pseudorapidity without observed hadrons) could be used for background sup- 
pression [31]. Unfortunately, in pp collisions of ^/s = 14 TeV at the LHC, overlapping 
events in a single bunch crossing will likely fill a rapidity gap even if it is present at the 
level of a single pp collision. Very low luminosity running is not useful because of the 
small signal cross section. 

The different color structures of signal and background processes can be exploited 
even at high luminosity, however, if one defines rapidity gaps in terms of minijets {ptj ~ 
15-40 GeV) instead of soft hadrons [32] . As has been shown for the analogous EW Zjj 
process [37], with its very similar kinematics, minijet emission in EW exchange occurs 
mainly in the very forward and very backward regions, and even here is substantially 
softer than in the QCD Zjj analogue. This is shown in Fig. 2.7, where we impose the 
cuts of Eqs. (2.1-2.3) and additionally require 

PTj > 40 GeV ; pre > 20 GeV , \r]e\ < 2 , (2.4) 

and 

- 10 GeV <m£^<M2 + 10 GeV. (2.5) 

Fig. 2.7 shows the rapidity distribution of the additional radiation (jet of pr > 20 GeV) 
with respect to the center of the two tagging jets. It is clear that the additional radiation 
in QCD processes tends to be much more central than in EW processes; the knee in the 
EW Zjj curve at \r]\ fti 2.5 is a phase space effect: that is the average rapidity of the 
tagging jets, and loss of rate is experienced due to the finite cone sizes of the tagging 
and third jets. This suggests that a veto on these central minijets will substantially 
improve the signal-to-background ratio. Following the analysis of Ref. [37], a veto is 
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Figure 2.7: Characteristics of the third (soft) jet in EW Zjjj "signal" (sohd hnes) and 
QCD Zjjj "background" (dashed hnes) events at the LHC. Acceptance cuts are as 
given in the text. The pseudorapidity rj^ is measured with respect to the center of the 



two tagging jets, rj = ^{rjf'^ + 77^"^ ) and the distributions are normahzed to the total 
respective cross sections. 



imposed on additional central jets, i.e., in the region 



PTi° > PT,veto ) 



tag veto ^ tag 

.max ' 



(2.6) 
(2.7) 



where PT,vcto may be chosen based on the capability of the detector; it is expected to 
be about 20 GeV for the LHC. 

Sizable background reduction via a minijet veto requires lowering of the pr,veto 
threshold to a range where the probability for additional parton emission becomes order 
unity. In a pcrturbative calculation the resulting condition, a{n + 1 jets) ^ a{n jets), 
indicates that one is leaving the validity range of fixed-order perturbation theory, and 
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it becomes difficult to provide refiable theoretical estimates of minijet emission rates. 
This happens at much higher values of pT,veto for QCD processes than for their EW 
analogues. Gluon emission is governed by very different scales in EW as compared to 
QCD processes, due to their different color structures. Thus, a parton shower approach 
does not immediately give reliable answers unless both color coherence and the choice 
of scale are implemented correctly, corresponding to the answer given by a complete 
QCD calculation. 

For the study of a central jet veto, the emission of at least one extra parton must 
be simulated. This is achieved by calculating the cross sections for the relevant "2-jet" 
process under consideration with additional gluon radiation and all crossing related 
processes. For the Hjj signal this includes 

qq qqHg , qq qqHg , qg qqqH , qg qqqH , (2.8) 

which can be found in Ref. [38]. For this case with three final-state partons, the factor- 
ization scale is chosen as /if — min(pT) of the tagging jets and the renormalization scale 
Hr is set to the transverse momentum of the non-tagging parton (minijet). Different 
scale choices or different input parameters will, of course, affect our numerical results. 

While the necessary information on angular distributions and hardness of additional 
radiation is available in the "3-jet" and tt+jets processes discussed in Chapters 3-5, one 
must either regulate or reinterpret these divergent cross sections. I use the truncated 
shower approximation (TSA) [39] for the former, treating the "2-jet" cross sections 
as the inclusive rate. For the latter I assume that additional soft parton emission, 
which will be dominated by soft gluons, exponentiates like soft photon emission. A 
Poisson distribution can then be used to estimate the veto probability. Details of these 
procedures can be found in Appendix C and Refs. [35,37]. Once the veto probability 
Pveto for a given process is known, the total rate is estimated by multiplying this by the 
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Table 2.2: Summary of veto survival probabilities for = 20 GeV used in Chap- 
ters 3-5. 



search 


Hjj 


tt 


tij, 


QCD 


EW 


QCD 


QCD 


DPS 








ttjj 


V{V)jj 


V{V)jj 


Wjjj 


bbjj 


lljj 




0.89 






0.30 


0.75 






0.30 




0.89 


0.46 


0.29 


0.29 


0.75 








rrjj 


0.87 






0.28 


0.80 


0.28 


0.28 





"2-jet" cross section, 

^surv — ^2— jet ' (1 Pveto) — ^2— jet ' Psurv • (2"^) 

The values for Pgurv I use in the following Chapters for the various processes considered 
are shown in Table 2.2. 

2.4 Summary 

I have identified and discussed the important characteristics of the core signal process, 
H production via weak boson fusion, in association with two quark jets, and estab- 
lished base kinematic cuts useful in enhancing the signal relative to a large class of 
backgrounds. These are the far forward and backward tagging jets of moderate pxi and 
a centrally produced Higgs with considerable pt, which will tend to yield decay prod- 
ucts with relatively high px- Additional cuts will be necessary to yield an observation 
of the Higgs, but will uniquely depend on the decay final state of the Higgs, and are 
discussed in the following Chapters. 

I have also identified an additional important distinction between the signal WBF 
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process and general QCD backgrounds, the collective hardness and rapidity distribu- 
tions of additional partonic activity in the central region between the tagging jets, 
attributable to higher-order gluon emission from the hard scattering process. I have 
discussed the method for calculating these distributions, which are only marginally per- 
turbative, mentioned two techniques for extracting a perturbative result, and provided 
a summary of the results. The techniques are discussed in detail in Appendix C. 
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Chapter 3 
The Search for ^ 77 

3.1 Introduction 

Fits to precision electroweak (EW) data have for some time suggested a relatively low 
Higgs boson mass, in the 100 GeV range [40] and this is one of the reasons why the 
search for an intermediate mass Higgs boson is particularly important [23]. Beyond 
the reach of LEP, for masses in the 100 — 150 GeV range, the H — > 77 decay channel 
at the CERN LHC is very promising. Consequently, LHC detectors are designed with 
excellent photon detection capabihties, resulting in the case of the CMS detector in a 
di-photon mass resolution of order 1 GeV for a Higgs boson mass around 120 GeV [41]. 
Another advantage of the i/ ^ 77 channel, in particular compared to the dominant 
H ^ bb mode, is the lower background from QCD processes. 

In this Chapter I demonstrate the observation potential of if —> 77 in WBF, out- 
lining the additional cuts needed to enhance the signal relative to the principal back- 
grounds, and including estimates of minijet veto survival probability. This discussion 
closely follows that in Ref. [34], but here I reanalyze this mode to include detector 
resolution effects, identification efficiencies and improved mass resolution estimates, 
corrections to the if — 77 branching ratio, and estimates of the minijet veto survival 
probabilities, found in Appendix C. A programming bug was also fixed that affected 
the QCD and EW 77jj backgrounds, increasing them by 20%. 
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3.2 Calculational Tools 

Higgs production and its characteristics has already been described in Section 2.2. We 
must now consider Higgs decay to the final state 77. As the Higgs is a spin-0 particle, 
it is sufiicient to treat its decay simply by multiplying the total cross section by the 
appropriate branching ratio and using phase space distributions for the decay products. 
For the H 77 partial decay width I use input from the program hdec AY [42] . In the 
mass region of interest, 100-150 GeV, B[H 77) varies from 0.13% to 0.22%, with 
the maximum value at about 125 GeV. As the photons are electromagnetic in nature, 
their momenta are smeared according to the resolution prescribed by Eq. B.2. 

In addition to the basic features of the Hjj signal, a centrally produced Higgs 
and two far forward, semi-hard tagging jets, the i7 — > 77 process will produce hard, 
central, isolated photons. These additional characteristics are shown in Fig. 3.1. The 
two photons tend to be emitted with considerable pr, in the 20-50 GeV range, due to 
the massive nature of the Higgs, relative to massless photons. 

Given the features of the signal, one must consider background processes that can 
lead to events with two hard, isolated photons and two forward jets. The largest 
background consists of all QCD 2 — 2 processes which contain one or two quark lines, 
from which the two photons are radiated. Examples are qQ ?Q77 or qg — > gg'77. 
Matrix elements for the irreducible QCD processes are available in the hterature [43]. 
For this irreducible QCD background, the renormalization scale and factorization scales 
are chosen as described in Appendix A. 

The EW analog to the above QCD processes is an additional background. The ir- 
reducible EW background consists of qQ — > qQ processes mediated by i-channel Z, 7, 
or W exchange, with additional radiation of two photons. 7 and Z exchange processes 
have amplitudes which are proportional to those of analogous gluon exchange processes, 
but with smaller couplings. These are ignored because, in all regions of phase space. 
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Figure 3.1: Normalized transverse momentum distributions of (a) the minimum-p^^ pho- 
ton and (b) the maximum-p^ photon in jj'j'j events at the LHC. The core acceptance 
cuts of Eqs. (2.1-2.3) are imposed. Results are shown for the qq — > qqH signal at 
niH — 120 GeV (sohd hne), the irreducible QCD background (dashed hne), the irre- 
ducible EW background (dot-dashed line), and for the double parton scattering (DPS) 
background (dotted line). 

they constitute only a tiny correction to the irreducible QCD background. All charged 
current qQ — > qQll (and crossing related) processes arc included, however; while for- 
mally of order and thus suppressed compared to the order Higgs signal, the small 
— > 77 branching ratio leads to comparable event rates. Furthermore, W exchange 
processes can produce central photons by emission from the exchanged W and, there- 
fore, are kinematically similar to the signal. This signal-like component remains after 
forward jet tagging cuts, as can readily be seen in the pt distribution of the jets in 
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Figure 3.2: Transverse momentum distributions of (a) the softer and (b) the harder of 
the two tagging jets in jj'yy events. Generic acceptance cuts (Eq. (2.1)) and forward 
jet tagging cuts (Eq. (2.2)) are imposed. 

Fig. 3.2. I am not aware of previous calculations of the irreducible EW background, 
and construct the matrix elements for it using madgraph [4]. The irreducible EW 
background is determined with the same choice of factorization scale as the signal. 

Double parton scattering (DPS), with pairs of jets and/or photons arising from two 
independent partonic collisions in one pp interaction, must also be considered. I do not, 
however, consider reducible backgrounds, where e.g. a jet fragmenting into a leading 7r° 
is misidentified as a photon. Reducible backgrounds were shown to be small compared 
to irreducible ones in the analysis of the i7 77 signal [41] and I assume the 

same to hold for the cleaner signal considered here. 

With jet transverse momenta as low as 20 GeV, double parton scattering (DPS) is 
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a potential source of backgrounds that must also be considered. DPS is the occurrence 
of two distinct hard scatterings in the collision of a single pair of protons. Following 
Ref. [44], I calculate the cross section for two distinguishable processes happening in 
one pp collision, as 

f^DPS = , I'J.ij 

with the additional constraint that the sum of initial parton energies from one proton 
be bounded by the beam energy, a^ff parameterizes the transverse size of the proton. 
It has been measured by the Fermilab CDF collaboration to be a^/f = 14.5 mb [45]. I 
assume the same value to hold for LHC energies. 

One DPS background arises from simultaneous 77J and jj events, where the jet in 
the 77J hard scattering is observed as a tagging jet, together with one of the two jets 
in the dijet process. In order to avoid a three-jet signature, one might want to require 
the second jet of the dijet process to fall outside the acceptance region of Eq. (2.1). 
However, this would severely underestimate this DPS contribution, since soft gluon 
radiation must be taken into account in a more realistic simulation. Soft radiation 
destroys the pr balance of the two jets in the dijet process, leading to the possibility of 
only one of the two final state partons to be identified as a jet, even though both satisfy 
the pseudo-rapidity requirements of Eq. (2.1). Since our tree- level calculation cannot 
properly take into account such effects, I conservatively select the higher-energy jet of 
the dijet process in the hemisphere opposite that of the jet from the 77J event, and 
allow the third jet to be anywhere, completely ignoring it for the purposes of imposing 
further cuts. I choose the factorization scale for DPS to be Uf = YIpt „■, and the 
renormalization scale to be Lir — J2pt ■ 

A second DPS mode consists of two overlapping events. All final-state particles 
must be observed above threshold in the detector. With full acceptance cuts this back- 
ground is found to be insignificant compared to the others, and will not be considered 
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further. I do not consider DPS backgrounds from overlapping 77 and jj events since 
the double forward jet tagging requirements of Eq. (2.3) forces the dijet invariant mass 
to be very large, effectively eliminating this background. 

3.3 Separation of Signal and Background 

As a reference starting point, our signal cross section times branching ratio with the 
cuts of Eqs. (2.1-2.3) for run = 120 GeV is anjj ■ B{H 77) = 2.4 fb. To ensure 
a clean signal the search must first establish observation of the Higgs decay products 
well-separated from the tagging jets. Thus our first additional cut, beyond the basic 
requirements of Eqs. (2.1-2.3), is a minimum px requirement for the photons: 



This still leaves a 2.2 fb signal for mu = 120 GeV, better than 90% acceptance. The 
backgrounds, however, are still overwhelming, as shown in the first line of Table 3.1. 
This is not surprising, as the presence of pt — 20 GeV jets is a common occurrence in 
hard scattering events at the LHC. 

Prominent features of the irreducible QCD background are the steeply falling trans- 
verse momentum distributions of both the jets and photons, as shown in Figs. 3.1 
and 3.2. These distributions are typical for bremsstrahlung processes and the low 
nature of the backgrounds is enhanced by the soft photons being emitted typically at 
high rapidity, forward of their parents quark jet. This allows one to suppress the back- 
grounds further by harder pt cuts; I find that the following asymmetric pt cuts bring 
the backgrounds down another factor of three, while accepting ^ 85% of the signal: 



Pt^ > 20 GeV . 



(3.2) 



Pt,, > 40 GeV , 



Pt,, > 20 GeV , 



PT,, > 50 GeV , PT,, > 25 GeV . 



(3.3) 
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For niH — 120 GeV, the resulting signal cross section is 1.9 fb. One could require even 
higher pt^ cuts, making the backgrounds negligible, but this would come at the expense 
of a sizeable reduction in signal rate. 

The SM Higgs resonance for this mass range is extremely narrow, 18 MeV at the 
most, which is much smaller than either of the expected photon resolutions of the CMS 
or ATLAS detectors, ~ ±0.6 ■ ■ ■ 1 GeV and ~ ±1.5 • • • 1.8 GeV, respectively. I use the 
expected resolution of the CMS detector, and examine 2 GeV mass bins: 

TTiH-l GeV < m^^ < TTiH + I GeV . (3.4) 

For this optimization, 70% of the Higgs signal's Gaussian profile would be retained [46]. 
This effect is included on line three of Table 3.1. Line four of the same table further 
introduces common CMS and ATLAS efficiencies for identification of tagging jets and 
photons: each tagging jet is a factor 0.86 and each photon is a factor 0.80, for both the 
Higgs signal and all backgrounds. 

Including these detector efficiencies and resolutions is an improvement on the VV —>■ 
if — > 77 analysis first presented in Ref. [34]. I also correct here a programming error 
that affected the QCD and EW 77JJ cross sections; the result is an approximate 20% 
increase in rate for those processes. The efficiencies and mass resolution demonstrate an 
overall signal rate loss of about two-thirds, which can fortunately be partly compensated 
for by exploiting another feature of the qq — > qqH signal, namely the absence of color 
exchange between the two scattering quarks. From the discussion of the minijet veto 
in Section 2.3 and summarized results for the Higgs signal, and general QCD and EW 
backgrounds using the TSA and detailed in Appendix C, I apply expected minijet veto 
survival probabilities to the signal and background rates, summarized in line five of 
Table 3.1. The veto is on additional radiation in the region between the two tagging 
jets, with pt > 20 GeV. 

The minijet veto improves signal-to-background (S/B) rates to better than 1/1 in all 
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Table 3.1: Signal tuh = 120 GeV and background ■yjjj cross sections (fb) for successive 
levels of cuts given by the Equations in parenthesis, apphcation of CMS expected tagging 
jet and photon ID efficiencies, mass resolution effects for 2 GeV bins, and apphcation 
of a minijet veto with = 20 GeV. 



cuts 




QCD jj77 


EW jj77 


DPS 


forward tagging + ID (2.1-2.3,3.2) 


2.2 


215 


62 


83 


+ staggered PT(j,'y) (3.3) 


1.9 


66 


29 


17 


+ 2 GeV mass bin (3.4) 


1.3 


0.87 


0.34 


0.24 


+ efficiencies (e = 0.473) 


0.63 


0.41 


0.16 


0.12 


Psurv,20 


xO.89 


xO.30 


xO.80 


xO.30 


+ minijet veto (2.6) 


0.56 


0.12 


0.12 


0.04 



cases, and better than 2/1 over much of the mass range examined. Gaussian equivalent 
statistical significances for Poisson statistical treatment of background fiuctuations is 
at the 5(7 level for almost the entire mass range with about 40-50 fb~^ of data at low 
luminosity. The shghtly lower significances at the lower end of the mass range are of 
no concern, as this region will already have been explored by both the CERN LEP and 
Fermilab Tevatron Higgs searches. 

3.4 Discussion 

Fig. 3.3 shows the results after the application of all cuts, Eqs. (2.1-2.3,3.2-3.4), inclusion 
of tagging jet and photon ID efficiencies, mass resolution effects for 2 GeV bins, and a 
minijet veto with pT,veto = 20 GeV. This plot compares the total signal cross section. 
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Table 3.2: Signal and total background 77jj cross sections (fb) for various Higgs masses, 
after application of all cuts, including a 2 GeV mass bin centered around the ex- 
pected Higgs mass, application of CMS expected tagging jet and photon ID efficiencies 
(e = (0.86)^ ■ (0.8)^ = 0.473), mass resolution effects for 2 GeV bins, and application 
of a minijet veto with p^*^" — 20 GeV. Gaussian equivalent Poisson statistical signal 
significances are given for 50 fb ~^ of data at low luminosity. 



Higgs mass (GoV) 


100 


110 


120 


130 


110 


150 


e-aHjj-B{H ^^^) (fb) 


0.37 


0.48 


0.56 


0.56 


0.48 


0.33 


e • ctqcd (fb) 


0.14 


0.13 


0.12 


0.11 


0.10 


0.08 


e • aEw (fb) 


0.14 


0.13 


0.12 


0.11 


0.10 


0.09 


e • aops (fb) 


0.05 


0.04 


0.04 


0.03 


0.03 


0.02 


e ■ (^bkg,tot (fb) 


0.33 


0.31 


0.28 


0.25 


0.22 


0.19 


S/B 


1.1 


1.6 


2.0 


2.3 


2.2 


1.8 


(T'Gaus 


3.8 


5.0 


6.0 


6.2 


5.7 


4.3 



in fb, to the di-photon invariant mass distribution, e • da/dm^^ in fb/GeV and thus 
indicates the relative size of signal and background. For our cuts, with 50 fb~^ of 
data at low luminosity, we thus expect anywhere from 16 to 28 i/ ^ 77 events on a 
background of 9 to 14 events, corresponding to a 4.3 to 6.2 standard deviation signal, 
depending on the Higgs mass. Thus, Higgs observation with a modest 40-50 fb~^ of data 
appears quite feasible in the qq qqH — •yyjj channel. We note that this channel 
is thus competitive with the inclusive i7 — > 77 search, which is predicted to require 
about 20 fb~^ (> 100 fb~^) of data to reach 5a significance coverage over the mass 
range 100-150 GeV by the CMS (ATLAS) collaboration [41]. 

A more detailed analysis is warranted because more than 50% of the signal events 
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have at least one jet with Irj] < 2.4 (see Fig. 2.3), leading to charged particle tracks 
in the central detector. As a result, the position of the interaction vertex can be more 
accurately obtained, leading to improved photon invariant mass resolution. We leave 
detailed studies of detector performance to the detector collaborations. 
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Figure 3.3: Higgs signal cross section (fb) and diphoton invariant mass distribution (in 
fb/GeV) for the backgrounds after application of all cuts, Eqs. (2.1-2.3,3.2, 3.3), includ- 
ing CMS expected tagging jet and photon ID efficiencies (e = (0.86)^ • (0.8)^ = 0.473), 
mass resolution effects, and application of a minijet veto with p^^° = 20 GeV. The 
squares are the Higgs signal for = 100, 110, 120, 130, 140, 150 GeV. The sohd hue 
represents the sum of all backgrounds, with individual components from the irreducible 
QCD background (dashed hne), the irreducible EW background (dot-dashed line), and 
the DPS background (dotted line) shown below. 
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Chapter 4 
The Search for H 

4.1 Introduction 

While I have estabhshed an additional technique to search for the Higgs in the mass 
range 100 < mn < 150 GeV, just above this range the branching fraction B{H 77) 
falls off sharply due to the increasing availabihty of the H W^*^W^*^ mode which 
dominates - the photon decay mode is no longer viable, even with very large luminosity. 
A search for the very clean four-lepton signature from H ZZ decay can find a Higgs 
boson in the mass region 130 ^ mn ^190 GeV, but due to the small branching fraction 
of this mode, very large integrated luminosities, up to 100 fb~^ or more, are required. 
Since H — > decays dominate in this region, I turn our attention to the cleanest 

WW decay signature, e^iJ^fr- One can search for this signature in inclusive gg H 
Higgs production, and while some attention has been given to this for the LHC [47,48] , 
production via WBF for the same decay mode has only recently been discussed by us 
in the hterature [36] . This discussion closely follows that work, but includes additional 
detector efficiencies for tagging jet identification. I will show that WBF is comparable 
to gluon fusion for discovery in this mode, and that it overlaps considerably with the 
photon decay region, allowing for additional coupling ratios to be determined with low 
luminosity, on the order of 10 fb~^. 
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In Section 4.2 I describe the calculational tools, the methods employed in the sim- 
ulation of the various processes, and important parameters. Extra minijet activity is 
simulated by adding the emission of one extra parton to the basic signal and background 
processes. Generically I call the basic signal process (with its two forward tagging jets) 
and the corresponding background calculations "2-jet" processes, and refer to the sim- 
ulations with one extra parton as "3-jet" processes. In Section 4.3, using the 2-jet 
programs for the backgrounds, I demonstrate forward jet tagging, a b veto and other 
important cuts which combine to yield an 2/1 to 1/2 signal-to-background (S/B) 
ratio, depending on the Higgs mass. 

In Section 4.4 I analyze the different minijet patterns in signal and background, using 
both the truncated shower approximation (TSA) [39] to regulate the cross sections, and 
the gluon exponentiation model to estimate the minijet multiplicity [57]. By exploiting 
the two most important characteristics of the extra radiation, its angular distribution 
and its hardness, the QCD backgrounds can be suppressed substantially by a veto on 
extra central jet emission. Within the TSA and exponentiation models, probabilities 
are estimated for vetoing signal and background events, and are combined with the 
production cross sections of the previous section to predict signal and background rates 
in Table 4.2. These rates demonstrate the feasibility of extracting a very low background 
H W^*^W^*^ signal at the LHC. 

The signal selection is not necessarily optimized yet. The variables I identify for cuts 
are the most distinctive, but deserve a multivariate analysis with detector simulation. 
I do demonstrate an additional variable in Section 4.5 which is not used for cuts, but 
rather can be used to extract the Higgs boson mass from the final event sample. 
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4.2 Calculational Tools 

All signal and background cross sections are determined in terms of full tree level matrix 
elements for the contributing subprocesses and are discussed in more detail below. 

For all our numerical results I have chosen the values listed in Appendix A, includ- 
ing the HDECAY-corrected effective branching ratios presented in Table A.l. Unless 
otherwise noted the factorization scale is chosen as fxj = min(pT) of the defined jets. 

4.2.1 The qq qqH(g),H W^*^W^*^ e^ii^i)T signal process 

An important additional tool for distinguishing the H — > e^is^^T signal from various 
backgrounds is the anti-correlation of the W spins, as pointed out in Ref. [48]. This is 
due to the preservation of angular momentum in the decay of the spin-0 Higgs boson. 
Of course, one can observe only the angular distributions of the charged decay leptons, 
but this is sufficient. The decay rate is proportional to (p^- ■ Pi^){pe+ ■ Pv)- In the 
rest frame of the Higgs boson, in which the e~P or e^v pairs are emitted back-to- 
back for W'^W' production at threshold, this product is a maximum for the charged 
leptons being emitted parallel. This characteristic is preserved and even enhanced when 
boosted to the lab frame, as the Higgs boson in weak boson fusion is typically emitted 
with pt ~ 60 - 120 GeV. 

To model this one cannot simply apply the H WW branching ratio to our produc- 
tion cross section and generate the Higgs decay products by phase space distributions 
as in the if — > 77 case. Instead, one must calculate the value of the matrix element 
squared for H WW e^iJ^-^T) separately for each phase space point. This re- 
sult can then be multiplied directly into the production cross section, thus maintaining 
angular correlations of the W decays. 
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4.2.2 The QCD + jets backgrounds 

Given the H decay signature, the main physics background to our e^//^^T signal arises 
from tt+jets production, due to the large top production cross section at the LHC and 
because the branching ratio B{t ^ Wb) is essentially 100%. 

The basic process I consider is pp tt, which can be either gg- or grg-initiated, with 
the former strongly dominating at the LHC. QCD corrections to this lead to additional 
real parton emission, i.e., to ti+j events. Relevant subprocesses are 

gq^ttq, gq ^ ttq , qq ^ ttg , gg ^ ttg , (4.1) 

and the subprocesses for tt + jj events can be obtained similarly. For the case of 
no additional partons, the 6's from the decaying top quarks may be identified as the 
tagging jets. In this case, calculating the cross section for tt + j where the 6's are 
explicitly identified as the tagging jets serves to estimate the effect of additional soft 
parton emission, i.e., minijct activity in the central detector; this is described in detail 
in Sec. 4.4. At the same time, one can identify a distinctly different, perturbative region 
of phase space, where the final-state light quark or gluon gives rise to one tagging jet, 
and one of the two decay 6's is identified as the other tagging jet. In this case, tt + jj 
may be used to estimate minijet activity for the hard process pp — > ti+j. Finally, there 
is a third distinct region of phase space, for the perturbative hard process pp tt + jj, 
where the final state light quarks or gluons are the two tagging jets. 

Thus, the "ttj" and "ttjj" calculations serve a dual purpose: to obtain the cross 
sections for the contribution of the perturbative processes where light quark or gluon 
jets lie in the region of phase space where they are experimentally identified as far- 
forward/backward tagging jets; and to estimate the additional QCD radiation patterns 
for the ncxt-lower-order perturbative tt + jets process. The tt and ttj matrix elements 
were constructed using madgraph [4] , while the ttjj matrix elements are from Ref . [49] . 
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Decays of the top quarks and VF's are included in the matrix elements; however, 
while the W's arc allowed to be off-shell, the top quarks are required to be on-shell. 
Energy loss from b iuX is included to generate more accurate distributions. In 
all cases, the factorization scale is chosen as Hf — mm{ET) of the massless partons/top 
quarks. The overall strong coupling constant factors are taken as (as)"' = Y[i=i c^siExJ, 
where the product runs over all light quarks, gluons and top quarks; i. e. the transverse 
momentum of each additional parton is taken as the relevant scale for its production, 
irrespective of the hardness of the underlying scattering event. This procedure guaran- 
tees that the same factors are used for the hard part of a ti+jets event, independent 
of the number of additional minijets, and at the same time the small scales relevant for 
soft-parton emission are implemented. 

4.2.3 The QCD WW +jj background 

The next obvious background arises from real-emission QCD corrections to W'^W~ 
production. For W~^W~jj events these background processes include [43] 

qg qgW^W- , qq qq'W+W- , (4.2) 

which are dominated by t-channel gluon exchange, and all crossing related processes, 
such as 

qq ^ ggW^W- , gg ^ qqW^W' . (4.3) 

I call these processes collectively the "QCD WWjj" background. I do not calculate 
cross sections for the corresponding VFl^-|-3-jet processes, but instead apply the minijet 
veto probabihty estimate for QCD Zjj(j), Z ^ tt processes, for the acceptance cuts 
of this search, found in Appendix C. The QCD rrjj and QCD WWjj backgrounds 
are quite similar kinematically, which justifies the use of the same veto probabilities for 
central jets. 
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The factorization scale is chosen as for the Higgs boson signal. The strong couphng 
constant factor is taken as (tt^)^ = as{pTi)c(s{pT2)j '^•e-i the transverse momentum of 
each additional parton is taken as the relevant scale for its production. Variation of 
the scales by a factor of 2 or | reveals scale uncertainties of 35%, however, which 
emphasizes the need for experimental input or NLO calculations. 

The WW background lacks the marked anti-correlation of W spins seen in the signal. 
As a result the momenta of the charged decay leptons will be more widely separated 
than in if ^ W^*^W^*^ events. 

4.2.4 The EW WW + jj background 

These backgrounds arise from W'^W~ bremsstrahlung in quark- ( ant i) quark scattering 
via t-channel electroweak boson exchange, with subsequent decay W^W~ —>■ i'^i~;^T- 

qq qqW+W- (4.4) 

Naively, this EW background may be thought of as suppressed compared to the anal- 
ogous QCD process in Eq. (4.2). However, as in the analogous case for 77 jj events 
discussed in Chapter 3, it includes electroweak boson fusion, VV — > W'^W~ via s- or 
t-channel 7/Z-exchange or via WW 4-point vertices, which has a momentum and 
color structure identical to the signal. Thus, it cannot easily be suppressed via cuts. 

The matrix elements for these processes were constructed using madgraph [4]. I 
include charged-current (CC) and neutral-current (NC) processes, but discard s-channel 
EW boson and t-channel quark exchange processes as their contribution was found to 
be ~ 1% only, while adding significantly to the CPU time needed for the calculation. 
In general, for the regions of phase space containing far-forward and -backward tagging 
jets, s-channel processes are severely suppressed. I refer collectively to these processes 
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as the "EW WWjj" background. Both W's are allowed to be off-shell, and all off- 
resonance graphs are included. In addition, the Higgs boson graphs must be included 
to make the calculation well-behaved at large H^-pair invariant masses. However, these 
graphs include our signal processes and might lead to double counting. Thus, I set 
m// to 60 GeV in the EW WWjj background to remove their contribution. A clean 
separation of the Higgs boson signal and the EW WWjj background is possible because 
interference effects between the two are neghgible for the Higgs boson mass range of 
interest. 

Again an estimate of additional gluon radiation patterns is needed. This was first 
done for EW processes in Ref. [50], but for different cuts on the hard process. Here I 
use the results for EW Zjj{j), Z ^ tt processes found in Appendix C. The EW rrjj 
and EW WWjj backgrounds are quite similar kinematically, which justifies the use of 
the same veto probabilities for central jets. 

4.2.5 The QCD r+r" background 

The leptonic decay of r's provides a source of electrons, muons and neutrinos which 
can be misidentified as W decays. Thus, I need to study real-emission QCD corrections 
to the Drell-Yan process qq (2', 7) — > r+r". For T^T~jj events these background 
processes include [51] 

(19 (19^^^~ : qq' qq'r^T' , (4.5) 

which are dominated by t-channel gluon exchange, and all crossing-related processes, 
such as 

QQ 99't^'t' : 99 ^ qqr'^r" . (4.6) 

All interference effects between virtual photon and Z-exchange are included. I call 
these processes collectively the "QCD rrjj" background. The cross sections for the 
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corresponding rr + 3-jet processes, which I need for our modehng of minijet activity 
in the QCD rrjj background, have been calculated in Refs. [52-54]. Similar to the 
treatment of the signal processes, I use a parton-level Monte-Carlo program based on 
the work of Ref. [53] to model the QCD rrjj and rrjjj backgrounds. 

The factorization scale is chosen as for the Higgs boson signal. With n = 2 and 
n = 3 colored partons in the final state, the overall strong-coupling constant factors 
are taken as {cks)"' — IliLi '^s(PrJ, i-G- the transverse momentum of each additional 
parton is taken as the relevant scale for its production, irrespective of the hardness of 
the underlying scattering event. This procedure guarantees that the same factors are 
used for the hard part of a Zjj event, independent of the number of additional minijets, 
and at the same time the small scales relevant for soft-gluon emission are implemented. 

The momentum distributions for the r decay products are generated as for the Higgs 
boson signal. Because of the (axial-) vector coupling of the virtual Z, 7 to r's, the pro- 
duced T"*" and T~ have the same chirality. This correlation of the r polarizations is taken 
into account by calculating individual helicity amplitudes and folding the correspond- 
ing cross sections with the appropriate r+ and decay distributions, i.e. the full r 
polarization information is retained in the energy distribution of the r decay products. 

The dual leptonic decays of the r's are simulated by multiplying the T^T~jj cross 
section by a branching ratio factor of (0.3518)^/2 and by implementing coUinear tau 
decays with helicity correlations as discussed in Appendix D. I also use the results from 
this study for our minijet emission estimates, summarized in Appendix C. 

4.2.6 The EW r+r" background 

These backgrounds arise from Z and 7 bremsstrahlung in quark-(anti)quark scattering 
via t-channel electroweak boson exchange, with subsequent decay Z, 7 ^ t'^t": 

qq' qq'r^T' (4.7) 
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Naively, this EW background may be thought of as suppressed compared to the anal- 
ogous QCD process in Eq. (4.5). However, the EW background includes electroweak 
boson fusion, VV t~^t~, either via i-channel r/i^-exchange or via s-channel j/Z- 
exchange, and the latter has a momentum and color structure which is identical to the 
signal and cannot easily be suppressed via cuts. 

I use the results of Ref. [55] for our calculation which ignore s-channel EW boson 
exchange contributing to qq production, and Pauli interference of identical quarks. 
When requiring a large rapidity separation between the two quark jets (tagging jets) the 
resulting large dijet invariant mass severely suppresses any s-channel processes which 
might give rise to the dijet pair, and the very different phase space regions of the two 
scattered quarks make Pauli interference effects small. All charged-current (CC) and 
neutral-current (NC) subprocesses are included. The CC process dominates over NC 
exchange, however, mainly because of the larger coupling of the quarks to the W as 
compared to photon and Z interactions. I will refer to these EW processes as the "EW 
TTjf background. 

The r decay distributions arc generated according to the prescription in Appendix D. 
Since the programs of Ref. [55] generate polarization averaged r+r" cross sections, I 
must assume unpolarized r's. However, as for the QCD rrjj background, the t'^t" pair 
arises from virtual vector boson decay, resulting in a t+ and t~ of the same chirality. 
This correlation of the r polarizations is taken into account. 

In order to determine the minijet activity in the EW Zjj background I need to eval- 
uate the 0{as) real parton emission corrections. The corresponding 0{a^as) diagrams 
for 

qq' qq'g T+r~ , (4.8) 

and all crossing related subprocesses, have first been calculated in Ref. [37]. Production 
of the T-pair via Z and 7 exchange is considered. The factorization and renormalization 
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scales are chosen to be the same as for the Hjj signal, as this is also a hard EW process. 

I have previously examined other scale choices for the rr backgrounds [37], and found 
small uncertainties (~ ±15%) for the EW component, while variations for the QCD 
component reach a factor 1.5. I thus expect the signal and EW rrjj background cross 
sections to be fairly well determined at leading order, while the much larger theoretical 
uncertainty for the QCD rrjj background emphasizes the need for experimental input. 

4.3 Separation of Signal and Background 

The qq qqH, H W^*^W^*^ ji^vv dual leptonic decay signal is characterized 

by two forward jets and the W decay leptons (e,/x). To begin, I impose the core 
acceptance requirements of forward tagging found in Eqs. (2.1-2.3) and additionally a 
minimum observability cut on the leptons, 

Pti > 20GeV. (4.9) 

This ensures observation and isolation of the two tagging jets and Higgs final-state decay 
products, as well as imposes severe restrictions on the bremsstrahlung nature of the 
backgrounds: the generally higher rapidity of the W^s in the QCD WWjj background 
as compared to the Higgs signal, for example, is due to weak boson bremsstrahlung 
occurring at small angles with respect to the parent quarks, producing 1^'s forward of 
the jets. Line 1 of Table 4.1 shows the effect of these cuts on the signal and backgrounds 
for a SM Higgs boson of mass tuh = 160 GeV. 

Somewhat surprisingly, the EW WWjj background reaches 68% of the QCD WWjj 
background already at this level. This can be explained by the contribution from W, Z, 7 
exchange and fusion processes which can produce central W pairs and arc therefore 
kinematically similar to the signal. This signal-like component remains after the forward 
jet tagging cuts. 
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As is readily seen from the first line of Table 4.1, the most worrisome background 
is W pairs from ti + jets production. Of the 1080 fb at the basic cuts level, 12 fb are 
from tt, 310 fb are from ttj, and the remaining 760 fb arise from ttjj production. The 
additional jets (corresponding to massless partons) are required to be identified as far 
forward tagging jets. The ttjj cross section is largest because the tt pair is not required 
to have as large an invariant mass as in the first two cases, where one or both 6's from 
the decay of the top quarks are required to be the tagging jets. 

For the events where one or both of the 6's are not identified as the tagging jets, 
they will most frequently lie between the two tagging jets, in the region where 1 search 
for the W decay leptons. Vetoing events with these additional b jets provides a powerful 
suppression tool to control the top background. I discard all events where a 6 or 6 jet 
with pt > 20 GeV is observed in the gap region between the tagging jets, 

PTi, > 20GeV , rjj^rnin <Vb< Vj,max ■ (4-10) 

This leads to a reduction of tij events by a factor 7 while ttjj events are suppressed by 
a factor 100. This results in cross sections of 43 and 7.6 fb, respectively, at the level of 
the forward tagging + lepton ID cuts of Eqs. (2.1-2.3,4.9), which are now comparable 
to the other individual backgrounds. Note that the much higher b veto probability for 
ttjj events results in a lower cross section than that for ttj events, an ordering which 
will remain even after final cuts have been imposed (see below). 

QCD processes at hadron colliders typically occur at smaller invariant masses than 
EW processes, due to the dominance of gluons at small Feynman x in the incoming pro- 
tons. I observe this behavior here, as shown in Fig. 4.1. The three tt+jets backgrounds 
have been combined for clarity, even though their individual distributions are slightly 
different. I can thus significantly reduce much of the QCD background by imposing a 
lower bound on the invariant mass of the tagging jets: 

mjj > 650 GeV . (4.11) 
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Mjj (GeV) 

Figure 4.1: Normalized invariant mass distribution of the two tagging jets for the signal 
(solid) and various backgrounds: tt + jets (dotted), QCD WWjj (long dashed), EW 
WWjj (short dashed) , QCD rrjj (long dash-dotted) and EW rrjj (short dash-dotted) . 
The cuts of Eqs. (2.1-2.3, 4.9-4.10) are imposed. 

Another significant difference is the angular distribution of the charged decay lep- 
tons, and //^, relative to each other. In the case of the Higgs signal, the W spins 
are anti-correlated, so the leptons are preferentially emitted in the same direction, close 
to each other. A significant fraction of the various backgrounds does not have anti- 
correlated W spins. These differences are demonstrated in Fig. 4.2, which shows the 
azimuthal (transverse plane) opening angle, polar (lab) opening angle, and separation 
in the lego plot. I exploit these features by establishing the following lepton-pair angular 
cuts: 

(/.e^<105°, cos^eM>0.2, ARe^, < 2.2 . (4.12) 
It should be noted that while these cuts appear to be very conservative, for higher Higgs 
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Figure 4.2: Normalized angular distributions of the charged leptons: azimuthal opening 
angle, lab opening angle, and separation in the lego plot. Results are shown for a Higgs 
boson mass of 160 GeV and 190 GeV (sohd lines) and for the various backgrounds as in 
Fig. 4.1. Lepton angular separation is clearly smaller for the rriii = 160 GeV scenario. 
The cuts of Eqs. (2.1-2.3, 4.9-4.10) are imposed. 

boson masses the 0e/i and Ai?e// distribution broadens out to higher values, overlapping 
the backgrounds more. For ttih ~ 180 — 200 GeV these cuts are roughly optimized and 
further tightening would require greater integrated luminosity for discovery at this upper 
end of the mass range. Because of the excellent signal-to-background ratio achieved 
below, I prefer to work with uniform acceptance cuts, instead of optimizing the cuts for 
specific Higgs boson mass regions. 

I also examine the distributions for lepton-pair invariant mass, me/i, and maximum 
lepton pt, as shown in Fig. 4.3 for the case uih = 160, 190 GeV. As is readily seen, the 
QCD backgrounds and EW WWjj background prefer significantly higher values for 
both observables. Thus, in addition to the angular variables, I find it useful to restrict 
the individual pt of the leptons, as well as the invariant mass of the pair: 



ruea < 110 GeV , pt^< 120 GeV . 



(4.13) 
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Figure 4.3: Normalized distributions of the dilepton invariant mass and maximum 
charged lepton momentum after the cuts of Eqs. (2.1-2.3, 4.9-4.10). Results are shown 

for a Higgs boson mass of 160 GeV and 190 GeV (solid lines) and for the various back- 
grounds as in Fig. 4.1. The niH — 160 GeV curve peaks at lower values of rueix and 

,max 

These are particularly effective against the top backgrounds, where the large top mass 
allows for very high-pT leptons far from the tagging jets, and against the EW WWjj 
background, where the leptons tend to be well-separated in the lego plot. Again, the 
cuts are set quite conservatively so as not to bias a lower Higgs boson mass. Results 
after cuts (4.11-4.13) are shown on the third hne of Table 4.1, for the case of a 160 GeV 
Higgs boson. 

At this level of cuts the combined QCD and EW rrjj backgrounds exceed all other 
individual backgrounds, contributing 50% of the total. One can take advantage of the 
fact that in these backgrounds, the Z or 7 is emitted with quite high pt, on the order 
of 100 GeV, which contributes to large r boosts and causes the r decay products to be 
nearly coUinear in the lab frame. Within the coUinear approximation, the r momenta 
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Table 4.1: Signal rates cr • B{H —>■ c^/x^^t) for ttih = 160 GeV and corresponding 
background cross sections, in fb. Results are given for various levels of cuts and are 
labeled by equation numbers discussed in the text. The expected tagging jet identifica- 
tion efficiency is shown on line 5. In the last line the minijet veto in included. Line six 
gives the survival probabilities for each process, with p^f*" = 20 GeV. 



cuts 




tt + jets 


QCD 


EW 


QCD 


EW 


S/B 








WWjj 


WWjj 


TTjj 


TTjj 




forward tagging (2.1-2.3,4.9) 


17.1 


1080 


4.4 


3.0 


15.8 


0.8 


wl/65 


+ b veto (4.10) 
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1/5.1 


+ Mjj, ang. cuts (4.11-4.13) 


11.8 


2.8 


0.54 


0.50 


3.6 


0.4 


1.5/1 


-|- real r rejection (4.14) 


11.4 


2.6 


0.50 


0.45 


0.6 


0.08 


2.7/1 


+ tag ID efficiency {xO.74) 


8.4 


1.9 


0.37 


0.33 


0.45 


0.06 


2.7/1 


Psurv ^20 


xO.89 


xO.29 


xO.29 


xO.75 


xO.29 


xO.15 




+ minijet veto (2.6) 


7.5 


0.56 


0.11 


0.25 


0.13 


0.05 


6.9/1 



can be reconstructed knowing the charged lepton momenta and the missing transverse 

momentum vector [35,66]. Labeling by Xr^^Xr^ the fraction of r energy each charged 

lepton takes with it in the r decay, 'pT,x,'i>T,y can be used to solve the two equations 

(x,y transverse directions) for the two unknowns x^-i j (see Appendix D). For real r 

decays, the vector must lie between the two leptons, and apart from finite detector 

resolution the reconstruction must yield < x-n 2 < 1- For the Hjj signal and other 

backgrounds, the collinear approximation is not valid because the W^s receive modest 

boosts in the lab only. In this case, the vector will rarely lie between the two leptons, 

and an attempt to reconstruct a r pair will result in x^-^ < or < for 95% of the 

events ^ . Additionally, the "r pair" invariant mass that is reconstructed does not peak 
^Conversely, requiring > 0, > largely eliminates WW backgrounds and promises clean 
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at mz, even when it is positive. One can therefore apply a highly efficient cut against 
the QCD and EW rrjj backgrounds by vetoing events where an attempt to reconstruct 
a T pair in the coUinear decay approximation results in two "real" r's near the Z pole: 

Xt-i, > , mz - 25 GeV < jtIt-t < mz + 25 GeV . (4-14) 

The results of this final cut are shown in line four of Table 4.1. The r backgrounds are 
virtually eliminated, while the signal and the other backgrounds each lose fa 5%. 

4.4 Minijet Veto 

If we are to veto central b jets to reduce the ti+jets background to a manageable level, 
I must take care to correctly estimate higher-order additional central partonic emission 
in the signal and backgrounds. Fortunately, due to the absence of color exchange be- 
tween the two scattering quarks in EW processes, which includes our Hjj signal, one 
expects soft gluon emission mainly in the very forward and very backward directions. 
However, for QCD processes, which are dominated by ^-channel color octet exchange, 
soft gluon radiation occurs mainly in the central detector. Thus, when I estimate addi- 
tional central radiation with pr > 20 GeV to match our b veto condition, I will reject 
QCD background events with much higher probability than the EW processes. Our 
b veto is then automatically also a minijet veto, a tool for QCD background suppres- 
sion which has been previously studied in great detail for Hjj production at hadron 
colliders [32,35,50]. 

Following the analysis of Ref. [37] for the analogous EW Zjj process which would 
be used to "calibrate" the tool at the LHC, I veto additional central jets in the cen- 
tral region between the two tagging jets according to Eq. (2.6), again using the cutoff 

PT.veto — 20 GeV as for the if — > 77 case in Chapter 3. 
isolation of -ff —> rr —> e^/j.^'^T [56] . 
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While the necessary information on angular distributions and hardness of additional 
radiation is available in the "3-jet" and ti + jets processes discussed in Section 4.2, 
one must either regulate or reinterpret these divergent cross sections. As discussed in 
Section 2.3, 1 use the TSA [39] for the former, and directly apply the results, summarized 
in Appendix C, to the QCD and EW cases for the WWjj and rrjj backgrounds. For 
the ti+jets backgrounds, it is simpler instead to reinterpret the divergent higher-order 
cross sections using the exponentiation model, also discussed in Appendix C. However, 
the estimated value Pguw — 0.12 for the tt + jets background is much lower than 
the values experienced via either the TSA method or the exponentiation method for 
other QCD processes, and this difference is not yet fully understood. Hence, I apply 
conservatively the same value of Pgurv — 0.29 for this background as for the other 
QCD backgrounds. The veto probabilities are summarized in line six of Table 4.1. I 
emphasize that while these probabilities arc estimates only, they can be independently 
determined at the LHC in processes like Zjj and Wjj production [37,55]. For a Higgs 
mass of 160 GeV we are left with a signal cross section of 7.5 fb compared to a total 
background of 1.09 fb. 

4.5 Discussion 

So far I have considered a single Higgs boson mass of 160 GeV only. Since I have largely 
avoided mass-specific cuts, I can immediately extend our results to a larger range of 
niH- The expected number of signal events for 115 GeV < ttih < 200 GeV and an 
integrated luminosity of 5 fb~^ are shown in Table 4.2. For the same luminosity, 5.5 
background events are expected. Thus, the signal-to-background rate, S/B, is better 
than 1/1 for ^ 125 < Mh < 200 GeV, almost the entire mass range. In the second 
row of Table 4.2 the Poisson probabihties for this background to fluctuate up to the 
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Table 4.2: Number of expected events for the Hjj signal, for 5 fb~^ integrated lumi- 
nosity and application of all efficiency factors and cuts including a minijet veto, but for 
a range of Higgs boson masses. The total background is 5.5 events. As a measure of 
the Poisson probability of the background to fluctuate up to the signal level, the second 
line gives croauss, the number of Gaussian equivalent standard deviations. 



rriH 


115 


120 


130 


140 


150 


160 


170 


180 


190 


200 


no. events 


1.6 


3.6 


8.8 


15.8 


24.0 


37.5 


36.2 


29.9 


20.8 


16.3 


^Gauss 


0.6 


1.2 


3.0 


5.0 


7.1 


10.0 


9.8 


8.4 


6.3 


5.1 



signal level are given, in terms of the equivalent Gaussian significances which can be 
expected in the experiment, on average. For the mass range Mh = 155 — 180 GeV, our 
technique is slightly superior to that of gluon fusion [48] . The gluon fusion and WBF 
modes together should then be able to measure coupling ratios quite accurately with 
moderate additional luminosity. 

These results show that it is possible to isolate a virtually background free qq — > 
qqH, H WW signal at the LHC, with sufficiently large counting rate to obtain a 5cr 
signal (or much better) with a mere 5 fb~^ of data for the mass range 140-200 GeV. 
Extending the observability region down to 130 GeV requires at most 15 fb~^. This 
nicely overlaps the regions of observability for H —>■ ■yy (100-150 GeV) [34] discussed in 
Chapter 3. To reach 120 GeV would require 65 fb~^ at low luminosity, and to reach 
115 GeV would require f^i 165 fb"^. 

As the H — > WW mode is likely to be the discovery channel for the mass range 
130 GeV < niH < 200 GeV, we wish to be able to reconstruct the Higgs boson mass. 
At threshold, the two (virtual) 1^'s are at rest in the Higgs boson center-of-mass frame. 
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resulting in = iriup, so we can calculate the transverse energy of both the charged 
lepton and invisible neutrino systems, 



Using these results for the transverse energies, we may compute a transverse mass of 
the dilepton-|^j^ system. 



At threshold this is exactly the Higgs boson transverse mass. Below threshold, the 
relation = m,^p is still an excellent approximation, while above threshold it begins 
to lose validity as the W bosons acquire a non-zero velocity in the Higgs boson rest 
frame. But even at rriH = 200 GeV this "pseudo" transverse mass remains extremely 
useful for mass reconstruction. I show the dramatic results in Fig. 4.4, for Higgs boson 
masses of 130, 160 and 190 GeV. Clearly visible is the Jacobian peak at Mt^^ — mn, in 
particular ior mn — 160 GeV. The combined backgrounds are added to the Higgs signal, 
and are shown after application of all cuts and detector efficiencies, as well as both the b 
and minijet vetoes discussed in the previous Sections. The very low background, in the 
absence of a Higgs signal, is also shown. The high purity of the signal is made possible 
because the weak boson fusion process, together with the H W'^W~ e^jf^-^T 
decay, provides a complex signal with a multitude of characteristics which distinguish 
it from the various backgrounds. 

An important point to note is that once a Higgs signal has been identified and its 
mass approximately determined, even further enhancement of the signal can be had by 
binning events in M^^^^. This would be especially useful for a low-mass Higgs, since 
Fig. 4.4 shows clearly that the peak of the distribution for background events occurs at 




(4.15) 




(4.16) 
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160 GeV. 
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Figure 4.4: Dilepton-^r transverse mass distributions expected for a Higgs of mass mn 
= 130, 160, and 190 GeV (solid) after the cuts of Eqs. (2.1-2.3, 4.9-4.13) and application 
of all detector efficiencies (e = (0.86)^ = 0.74) and a minijet veto withpT^^veto — 20 GeV. 
Also shown is the background only (dashed). 

For H WW decays, lepton angular distributions are extremely useful for reducing 
the backgrounds even further. The anti-correlation of W spins in H decay forces the 
charged leptons to be preferentially emitted in the same direction, close together in the 
lego plot. This happens for a small fraction of the background only. I have identified 
the most important distributions for enhancing the signal relative to the background, 
and set the various cuts conservatively to avoid bias for a certain Higgs boson mass 
range. There is ample room for improvement of our results via a multivariate analysis 
of a complete set of signal and background distributions, which I encourage the LHC 
collaborations to pursue. Additional suppression of the tt + jets background may be 
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possible with b identification and veto in the pt < 20 GeV region. 

In addition to various invariant mass and angular cuts, 1 can differentiate between 
the Vl^'s of the signal and W, t backgrounds and the real r's in the QCD and EW 
TTjj backgrounds. This is possible because the high energy of the produced r's makes 
their decay products almost coUinear. Combined with the substantial px of the t'^t~ 
system this allows for r-pair mass reconstruction. The W decays do not exhibit this 
coUinearity due to their large mass, thus the angular correlation between the |ir vector 
and the charged lepton momenta is markedly different. Our real-r rejection makes use 
of these differences and promises to virtually eliminate the rrjj backgrounds. 

A final step is to veto additional soft jet activity in the central region between the 
two tagging jets, as additional gluon radiation in QCD processes is characteristically 
central and hard, whereas for WBF/EW processes it tends to be more soft and for- 
ward/backward. I expect a typical 70% reduction in QCD backgrounds for a central jet 
veto implemented for pt > 20 GeV, and about a 25% reduction for EW backgrounds 
but only about a 10% suppression for the WBF Higgs production. 
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Chapter 5 
The Search for H t^t~ 

5.1 Introduction 

For Higgs searches in the portion of the intermediate- mass range where decay to fermion 
pairs has a non-ncgUgible branching fraction, ^ 110 — 150 GeV, I show that observation 
of the H ^ TT decay channel is quite promising. An advantage of the H ^ tt channel, 
in particular compared to the dominant if — > 66 mode, is the lower background from 
QCD processes. The H ^ tt channel thus also offers the best prospects for a direct 
measurement of the Higgs boson's couplings to fermions. 

While some attention has been given to A/ H ^ rr searches at the LHC [41,59,60] 
in the framework of the MSSM, where the increased couphngs of A/H to r predicted 
for tan/9 ^ 1 lead to higher production rates, conventional wisdom says that the 
chance of seeing the SM Higgs via this decay mode is nil, and it had been ignored in 
the literature until our recent analysis [35] . I discuss those results here, and update the 
analysis with corrected H ^ rr branching ratios, inclusion of a tagging jet identification 
efficiency, and an improved estimate of minijet veto probabilities. The results apply to 
an intermediate-mass SM if —> rr at the LHC, covering the main physics and reducible 
backgrounds, and demonstrate the feasibility of Higgs boson detection in this channel 
with only modest luminosity, if — > rr event characteristics are analyzed for one r 
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decaying leptonically and the other decaying hadronically, because of the high trigger 
efficiency and good branching ratio of this mode; Ref. [59] found the dual leptonic decay 
mode to be considerably more difficult due to higher backgrounds. 

In Section 5.2 I describe our calculational tools, the methods employed in the sim- 
ulation of the various processes, and important parameters. Extra minijet activity is 
again simulated by adding the emission of one extra parton to the basic signal and 
background processes. In Sections 5.3 & 5.4, using the 2-jet programs for physics and 
reducible backgrounds, respectively, I demonstrate forward jet tagging and r identifica- 
tion and reconstruction criteria which yield an ?5i2/l signal-to-background (S/B) ratio. 
Both the Wj + jj and bbjj reducible backgrounds intrinsically are much larger than 
the Z — > TT and Drell-Yan r-pair production backgrounds. I explain and emphasize 
the cuts crucial to reducing these backgrounds to a manageable level. 

In Section 5.5 I analyze the different minijet patterns in signal and background, using 
the TSA [39] to regulate the cross sections. Within the TSA, probabilities are estimated 
for vetoing signal and background events, and are combined with the production cross 
sections of the previous section to predict signal and background rates in Table 5.2. 
These rates demonstrate the possibility to extract a very low background — > rr 
signal at the LHC. 

The signal selection is not necessarily optimized yet. Additional observables are 
available to distinguish the signal from background. The final discussion in Section 5.6 
includes a survey of distributions which can be used, e.g. in neural-net algorithms, to 
further improve the signal significance. 
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5.2 Calculational Tools 

Physical constants are chosen as in Appendix A, including HDECAY-corrected H ^ tt 
branching ratios. I employ CTEQ4L parton distribution functions [65] throughout. 
Unless otherwise noted the factorization scale is chosen as /if — min(pr) of the defined 
jets. 

5.2.1 The qq qqH{g); H t^t^ signal process 

In the following I consider only r-pair decays with one r decaying leptonically, r — > 
eUeUr, iii^hI't, and the other decaying hadronically, — > h^X, with a combined branch- 
ing fraction of 45%. 

Positive identification of the hadronic — > h^X decay requires severe cuts on the 
charged hadron isolation. Possible strategics have been analyzed by Cavalli et al. [59] 
and I base my simulations on their result. Considering hadronic jets of Et > 40 GeV 
in the ATLAS detector, they find non-tau rejection factors of 400 or more (see below) 
while retaining true hadronic r decays with an identification efficiency 

e^(T ^ u + hadrons) = 0.26 . (5.1) 

This estimate includes the requirement of seeing a single charged hadron track, of 
Pt > 2 GeV, pointing in the r direction, and thus effectively singles out 1-prong r 
decays. Accordingly, only the 1-prong hadronic branching ratios are considered in our 
mixture of tt, p and ai modes. Since the overall efficiency includes 3-prong events, which 
have negligible acceptance, the effective efficiency for 1-prong events is larger and taken 
as 0.337 in the following, which reproduces the overall efficiency of Eq. (5.1). 
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5.2.2 The QCD and EW r+r + physics backgrounds 

These backgrounds are identical to the QCD and EW rTjj{j) backgrounds discussed 
in Chapter 4, except that here I have one tau decaying leptonically and the other 
hadronically; the decays are discussed in Appendix D. 

5.2.3 The QCD Wj reducible background 

Reducible backgrounds to r+r~ events can arise from any process with a 

hard, isolated lepton, missing pr, and an additional narrow jet in the final state which 
can be mistaken as a hadronically decaying r. A primary reducible background thus 
arises from leptonic W decays in Wj events, where additional QCD radiation supplies 
the two tagging jet candidates. At lowest order I need to consider Wj+jj production as 
the hard process, which is very similar to the simulation of the QCD Zjjj background 
discussed before, with the bremsstrahlung Z replaced by a W. W — > ez/g, Z^^'^j decays 
only are considered and are treated as a fake r decaying leptonically. Real leptonic r 
decays from W — > Ti/r iveVrJ^r are relatively suppressed by the r leptonic branching 
ratio of 35% and the severity of the transverse momentum cuts on the softer charged 
lepton spectrum. They will be ignored in the following. 

Two of the jets in Wj + jj events are identified as tagging jets, and fluctuations 
of the third into a narrow jet are considered, resembling a hadronically-decaying r. In 
Ref. [59] the probability for misidentifying a gluon or light-quark jet as a hadronic r 
decay was estimated as 

e^(jet ^ "u + hadrons") = 0.0025 , (5.2) 

and I assign this probability to each of the final state jets. In each event one of the hard 
partons is randomly assigned to be the r. To mimic the signal, this jet and the identified 
charged lepton must be of opposite charge. Thus, I reduce the Wj + jj background by 
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an additional factor of two to simulate the opposite charge requirement for the single 
track allowed in the r-like jet. As the Wj + jj events are a QCD background, 1 use the 
same factorization and renormalization scales as for the QCD Zjj case. 

To simulate additional minijet emission, I need to add one more parton to the final 
state. The code for W+Aj matrix elements has been available since the work of Berends 
et al. [61]. Here I use the program developed in Ref. [62], which was generated via 
MADGRAPH [4] . Since W + Aj production produces a six-particle final state, with up to 
516 graphs for the most complicated processes, it takes considerable CPU time to obtain 
good statistics. I modified the madgraph code to do random helicity summation, 
speeding up the calculation by approximately a factor of 3 for a given statistical error 
in the final cross section. As before, ag is taken as the geometric mean of agipr) factors 
for each of the partons, including the parton which fakes the hadronic r decay. 

5.2.4 The QCD bbjj reducible background 

The semileptonic decay of 6-quarks provides another source of Icptons and neutrinos 
which can be misidentified as tau decays. Even though 6-quark decays are unlikely to 
lead to isolated charged leptons and very narrow tau-like jets in a single event, the sheer 
number of bb pairs produced at the LHC makes them potentially dangerous. Indeed, 
the analysis of Ref. [59] found that bb pairs lead to a reducible t~^t~ background which 
is similar in size to Wj production. I therefore study bbjj production as our second 
reducible background and neglect any other sources like tt events which were shown to 
give substantially smaller backgrounds to r^r^ -pairs in Ref. [59]. 

I only consider 6-production events where both 6-quarks have large transverse mo- 
mentum. In addition, two forward tagging jets will be required as part of the signal 
event selection. The relevant leading order process therefore is the production of bb 
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pairs in association with two jets, which includes the subprocesses 

gg bbgg 

qg bbqg (5.3) 
qiq2 bbqiq2 ■ 

The exact matrix elements for the 0{aj) processes arc evaluated, including all the 
crossing related subprocesses, and retaining a finite 6-quark mass [49] . The Pauli inter- 
ference terms between identical quark flavors in the process qiq2 bbqiq2 are neglected, 
with little effect in the overall cross section, due to the large differences in the rapidity 
of the final state partons. The factorization scale is chosen as the smallest transverse 
energy of the final state partons before the 6-quark decay. The strong coupling constant 
as is evaluated at the corresponding transverse energy of the final state partons, i.e., 
= as{ET{b))as{ETib))asipT,ict^)oisiPT,jetJ- 

The semileptonic decay b —>■ ivc of one of the 6-quarks is simulated by multiplying 
the bbjj cross section by a branching ratio factor of 0.395 (corresponding to at least 
one semileptonic 6-decay to occur) and by implementing a three-body phase space dis- 
tribution for the decay momenta. This part of the simulation is performed in order to 
estimate the effects of the lepton isolation cuts on the transverse momentum distribu- 
tions of the 6-decay leptons. Since these are kinematic effects I use the lightest meson 
masses in the simulation and set nib = 5.28 GeV and rric = 1.87 GeV. In Ref. [59] 
a factor 100 reduction of the bb background was found as a result of lepton isolation, 
requiring Et < 5 GeV in a cone of radius 0.6 around the charged lepton. In the sim- 
ulation, after energy smearing of the charm quark jet (see below), I find a reduction 
factor of 52 due to lepton isolation with a cone of radius 0.7. However, the simulation 
does not include parton showers or hadronization of the 6-quark, effectively replacing 
the 6-quark fragmentation function by a delta-function at one, and thus underestimates 
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the effect of lepton isolation cuts on the 6-quark background. Since I cannot model 
the lepton isolation efficiency, 1 multiply the bbjj rates by another factor 0.52, thus 
effectively implementing the factor 100 suppression found by Cavalli et al. [59]. This 
value should probably be reanalyzed in light of the different signal signature here. 

In addition to an isolated lepton, the bbjj events must produce a narrow jet which 
is consistent with a hadronic r decay, and has charge opposite the identified charged 
lepton. This may either be one of the light quark or gluon jets, for which the misiden- 
tification probability of 0.25% of Eq. (5.2) will be used, or it may be the 6-quark jet. 
In Ref. [59] the probability for misidentifying a 6-quark jet as a hadronic r decay was 
estimated as 

e^(6 ^ "v + hadrons") ^ 0.0005 . (5.4) 

However, due to limited Monte Carlo statistics, this number was based on a single 
surviving event only. Since we are really interested in an upper bound on the bbjj 
background I follow the ATLAS proposal [41] instead, and use the upper bound, 

er{b "v + hadrons") < 0.0015 , (5.5) 

for our analysis. Thus, all our bbjj cross sections, after r identification, should be 
considered conservative estimates. A more precise analysis of 6 ^ r misidentification 
probabilities in the LHC detectors is clearly needed, which is beyond the scope of the 
present work. Finally, an additional overall factor of two reduction is applied, as in the 
Wj + jj case, for the lepton-jet opposite charge requirement. 

The purpose of our 6-analysis is to verify that b semileptonic decays do not over- 
whelm the signal. The above procedures are adequate for this purpose, since I obtain 
final bbjj backgrounds (in Table 5.2) which are 20 to 40 times smaller than the signal. 
I do not calculate additional b quark backgrounds arising from intrinsic b contributions 
(processes like gb —>■ bggg). The matrix elements for these processes are of the same 
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order (a^) as for the bbjj subprocesses discussed above, but they are suppressed in ad- 
dition by the small 6-quark density in the proton. Also, I do not simulate additional soft 
gluon emission for the bbjj background. This would require 66 + 3 jet matrix elements 
which are not yet available. Rather, I assume the probability for extra minijet emission 
to be the same as for the other reducible QCD background, Wj + jj production. 

5.3 Higgs Signal and Real r+r~ Backgrounds 

I begin as in Chapters 3 & 4, adopting the core forward tagging requirements of 
Eqs. (2.1-2.3). As in the if — > 77 analysis I adopt asymmetric Pt cuts for the jets, 
a useful tool to discriminate against the steeply falling distributions typical of QCD 
backgrounds: 

> 40, 20 GeV . (5.6) 

The resulting Hjj, H ^ tt cross section is compared with the irreducible Zjj, Z — > 
TT backgrounds in the first row of Table 5.1. Somewhat surprisingly, the EW Zjj 
background reaches 5% of the QCD Zjj background already at this level, while naively 
one might expect suppression by a factor {olqed/ OLsf ~ 4 x 10~^. In the EW Zjj 
background, W exchange processes can produce central r pairs by Z emission from 
the exchanged W and are therefore kinematically similar to the signal. This signal-like 
component remains after the forward jet tagging cuts, and, as we will see, will grow in 
relative importance as the overall signal/background ratio is improved. 

So far I have not considered r decays. In order to get more realistic rate estimates 
and to include the reducible backgrounds {Wj + jj and bbjj, see Section 5.4) I need to 
study definite r decay channels. I consider r+r" decays with one r decaying leptoni- 
cally (e or //) and the other decaying hadronically in the following, since previous studies 
have shown that dual leptonic decay is more difficult to observe [59] . With a hadronic 
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Table 5.1: Signal and background cross sections Ba (fb) for mu = 120 GeV Hjj events. 
Results are given after increasingly stringent cuts given by the Equation numbers in 
parenthesis, and all values include the efficiency for tagging jet identification e = 0.74. 
The last column gives the ratio of the signal to the background cross sections listed in 
the previous columns. 



cuts 


Hjj 


QCD Zn 


EW Zjj 




bb + jj 


S/B 


forward tagging (2.1-2.3,5.6) 


50.6 


1240 


67 








+ r identification (5.8) 


1.47 


14.8 


1.07 


19.5 


5.6 


1/28 


+ 110 < rrirr < 130GeV (5.9) 


0.97 


0.70 


0.05 


1.31 


0.44 


1/2.6 


+ mjj > 1 TeV, mT{l,i)T) < 30 GeV 














(5.10,5.11) 


0.51 


0.12 


0.03 


0.08 


0.11 


1.5/1 


+ Xri < 0.75, Xr^ < 1.0 (5.12) 


0.40 


0.11 


0.02 


0.02 


0.04 


2.1/1 



branching ratio B{t —>■ u + hadrons) = 0.65 and the overall hadronic r-decay identi- 
fication efficiency of Eq. (5.1), the selection of this r-pair decay channel immediately 
reduces all t+t~ rates by a factor 

eB — 2€r{T u + hadrons) B{t ^ u + hadrons) B{t ^ ii'ei^r) 

= 2 •0.26- 0.65 -0.35 = 1/8.5. (5.7) 

In addition, triggering the event via the isolated r-decay lepton and identifying the 
hadronic r decay as discussed in Ref. [59] requires sizable transverse momenta for the 
observable r decay products. In the following I require 

PT.,., > 20 GeV , pT.,,„, > 40 GeV , (5.8) 

where the second requirement is needed to use the results of Ref. [59] on hadronic r 
identification. These transverse momentum requirements are quite severe and reduce 
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the Higgs signal by another factor of 3.8. Resulting signal and background cross sections 
are given in the second row of Table 5.1. 

Crucial for further background reduction is the observation that the r-pair invariant 
mass can be reconstructed from the observable r decay products and the missing trans- 
verse momentum vector of the event [63]. Details of the reconstruction are found in 
Appendix D. The assumption of coUinear tau decays is satisfied to an excellent degree 
because of the high r transverse momenta needed to satisfy Eq. (5.8). r pair mass recon- 
struction is possible only as long as the the decay products are not back-to-back. This 
last condition is met in our case because the H and Z bosons are typically produced 
with high on the order of 150 GeV for all processes except the bbjj background (in 
which case the 85 GeV is still sufficient). 

Mismeasured transverse momenta (smearing effects) can still lead to unphysical 
solutions for the reconstructed r momenta. In order to avoid these, 1 impose a cut on 
the angle between the r decay products and require positivity of the calculated x^-.: 

cos6't-t- > -0.9 , Xt-, ^>0. (5.9) 

The resulting r-pair invariant mass resolution is somewhat narrower than the one 
found in Ref. [59], the 1-a half- width for the H peak ranging from about 7 GeV for 
niH — 110 GeV to about 10 GeV for mn — 150 GeV (see Fig. 5.4 below). This improved 
resolution is an effect of the higher average pt of the underlying process: in our case, 
the two forward tagging jets from weak boson scattering impart a higher px on the H 
or Z than is the case from QCD radiation in gluon fusion. The smaller t^t~ opening 
angle then leads to a better r momentum reconstruction via Eq. (D.4). Given this 
r-pair mass resolution, I choose ±10 GeV mass bins for analyzing the cross sections. 
Signal and background cross sections in a 20 GeV mass bin centered at 120 GeV, after 
the reconstruction conditions of Eq. (5.9), are listed in the third row of Table 5.1. QCD 
and EW Zjj backgrounds are reduced by a factor of 20, while about 2/3 of the signal 
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Figure 5.1: Invariant mass distribution of the two tagging jets for the = 120 GeV 
Hjj signal (sohd hne) and the QCD Zjj background (dashed hne), at the level of 
forward tagging cuts and r reconstruction, Eqs. (2.1-2.3, 5.6-5.9,5.11). 

survives the mass reconstruction cuts. 

Because the QCD backgrounds typically occur at small invariant masses, I can 
further reduce them by imposing a cut on the invariant mass of the tagging jets. 



rrijj > 1 TeV. 



(5.10) 



Fig. 5.1 shows the tagging jets' invariant mass distribution for the signal and QCD Zjj 
background to illustrate the effect of the cut. 



5.4 Fake r+r Events: Reducible Backgrounds 

Reducible backgrounds to the H tt signal, with subsequent leptonic decay of one of 
the r's, arise from any source of isolated, single hard leptons. As discussed in Section 5.2, 
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I consider Wj + jj events and heavy quark production, in the form of bbjj events. 
Intrinsically, these reducible backgrounds are enormous and overwhelm even the physics 
backgrounds before r identification and tight lepton isolation cuts are made. Crucial 
for the reduction of these backgrounds to a manageable level is the requirement of a 
narrow r-like jet, which leads to a factor 400 suppression for the Wj+jj background (see 
Section 5.2.3). The probability for a 6-quark to fluctuate into a narrow r-like jet is even 
smaller, below 0.0015, and another large reduction, by a factor 100 (see Section 5.2.4), 
is expected from requiring the 6-decay lepton to be well isolated. An additional factor 
of two reduction is achieved by requiring opposite charges for the isolated lepton and 
the tau-like jet. The resulting background rates, for charged leptons and r-like jets 
satisfying the transverse momentum requirements of Eq. (5.8), are listed in the second 
row of Table 5.1. 

Unlike the Higgs signal or the Zjj backgrounds, the reducible backgrounds show 
no resonance peaks in the m^,- distribution. As a result, another reduction by an order 
of magnitude is achieved when comparing rates in a Higgs search bin of width 20 GeV 
(third row of Table 5.1). Additional reductions are possible by making use of specific 
properties of the reducible backgrounds. Analogous to the QCD Zjj background, the 
Wj + jj and bbjj backgrounds are created at smaller parton center of mass energies 
than the signal. As a result, the mjj > 1 TeV cut of Eq. (5.10) reduces both of them 
by roughly a factor of 4. 

Further suppression of the Wj +jj background can be achieved by taking advantage 
of the Jacobian peak in the lepton-i^T transverse mass distribution [59] , a feature which 
is otherwise used to measure the mass of the W. I compare the ttit distribution for the 
signal and the Wj + jj background in Fig. 5.2. A cut 

rriTiijT) < 30 GeV (5.11) 
reduces the Wj + jj background by a factor of 5 while reducing the signal acceptance 
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Figure 5.2: Transverse mass distribution of the 1-]/)t system for the Mh — 120 GeV 
Hjj signal (sohd hne) and the Wj+jj reducible background (dashed hne), at the level 
of far forward tagging cuts, r-reconstruction, and rrijj > 1 TeV (Eqs. 2.1-2.3, 5.6-5.10). 

by only 15%. Similar to the signal, the other backgrounds are affected very httle by 
the transverse mass cut. 

At this level the S/B ratio is nearly 1/1, and I study additional event characteristics, 
such as the missing momentum. As discussed previously in Section 4.3, in real r-pair 
events, the missing momentum is a vector combination of neutrino momenta, which 
carry away a significant fraction of the r"*" and r~ energies. In the reducible backgrounds 
it is purely from the leptonically decaying parent particle, either the W or one of the 6's. 
As such, one should reconstruct x^-^ = 1 for the narrow r-like jet, except for smearing 
effects. The effect is clearly observable in the distribution of events in the x^-Xr,^ plane, 
which is shown in Fig. 5.3. The Xn distribution of the leptonically decaying r-candidate 
also is softer for real r's than for the reducible backgrounds, because the charged lepton 
shares the parent r energy with two neutrinos. A cut 

Xri < 0.75 , Xr^ <1, (5.12) 
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proves very effective in suppressing tfie reducible backgrounds. For tfie Wj + jj back- 
ground 1 find suppression by another factor of 4.5 and tlie bbjj background is reduced 
by a factor of 3, wliile retaining 80% of tlie signal rate. One should note that these 
cuts are not optimized, they are merely chosen to demonstrate the usefulness of the 
Xti-Xt^ distributions in restricting the otherwise troublesome reducible backgrounds to 
a manageable level. Cross sections including these cuts are given in the last row of 
Table 5.1. 

In principle, the x-r distributions contain information on r polarization and Xti-Xt-^^ 
correlations allow one to distinguish between the decay of a spin-0 object, like the Higgs 
which results in opposite r"*" and r~ chiralities, and the decay of the spin-1 Z boson, with 
equal chiralities [64]. Comparison of the two scatter plots in Figs. 5.3(a) and 5.3(b) 
shows, however, that the remaining correlations are very weak. This may partially be 
due to the stringent transverse momentum cuts (5.8) on the r decay products which 
needed to be imposed for background reduction. In addition, the visible r energy 
fractions in r — > IVe^u-r and r — > pUr decays are mediocre polarization analyzers only 
(measuring the splitting of the p's energy between its two decay pions would improve 
the situation for the latter [58]). A dedicated study is needed to decide whether a r 
polarization analysis is feasible at the LHC, but because of the small rates implied by 
Table 5.1 I do not pursue this issue here. 

5.5 Minijet Veto 

Finally, I apply the results of a veto on additional central radiation, as for the if — > 77 
and H WW cases. Details of our calculations for the H ^ tt case may be found in 
Appendix C. In general, the QCD backgrounds here are rejected about three times as 
often as the Higgs signal, while the EW background is rejected only marginally more 



78 



1.50 
1.25 
1.00 
0.75 
0.50 
0.25 
0.00 
1.25 
1.00 
0.75 
0.50 
0.25 



0.00 



I ' ' ■ ■ I ' ■ ' ' i; ' ' ■ I 



-'•.4.1.. 1 I 



(a) 




(c) 




(d) 



0.00 0.25 0.50 0.75 1.00 1.25 0.00 0.25 0.50 0.75 1.00 1.25 1.50 



X 



Figure 5.3: Scatter plots of v. x^-^ with the cuts of Eqs. (2. 1-2. 3, 5. 8, 5.9-5.11), for: (a) 
the 120 GeV Hjj signal; (b) the combined QCD and EW Zjj irreducible backgrounds; 
(c) the Wj + jj and (d) the bbjj reducible backgrounds. The number of points in each 
plot is arbitrary and corresponds to significantly higher integrated luminosities than 
expected for the LHC. The solid hues indicate the cuts of Eq. (5.12). 

often, reflecting the component that shares the kinematic and color structure of the 
signal. 

Table 5.2 applies the survival probabilities found for the ?7-method of selecting tag- 
ging jets to the cross sections after final cuts, for Higgs boson masses ranging from 
110 to 150 GeV. A constant size of the mass bins of 20 GeV is kept for simplicity. 
In the actual experiment, the mass window will need to be optimized depending on 
the predicted width of the signal and background distributions, and may have to be 
asymmetric for low values of uih- Our table merely shows how observing a light Higgs 
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Table 5.2: Number of expected events for the signal and backgrounds, for 60 fb~^ 
at low luminosity and cuts as in the last line of Table 5.1 and a minijet veto with 
p>^*° — 20 GeV, including an efficiency factor for tagging jet identification (e — 0.74), 
for a range of Higgs boson masses. Mass bins of ±10 GeV around a given central value 
are assumed. 



mH(GeV) 


HjJ 


QCD Zjj 


EW Zjj 


Wj + jj 


bbjj 


^Gauss 


110 


24.2 


6.3 


3.4 


0.3 


0.8 


5.7 


120 


20.6 


1.8 


1.2 


0.3 


0.7 


7.4 


130 


16.0 


0.9 


0.7 


0.3 


0.6 


6.3 


140 


10.0 


0.6 


0.5 


0.4 


0.5 


4.7 


150 


4.8 


0.4 


0.4 


0.3 


0.4 


2.6 



boson is quite feasible, even in the mass window close to the smeared Z peak. As 
ruH approaches 150 GeV, however, the H ^ tt branching ratio drops rapidly in the 
SM and the signal gets low for integrated luminosities of order 60 fb~^ at low machine 
luminosity. It should be noted that with higher integrated luminosity, this channel is 
still very effective to make a direct measurement of the Htt coupling; it would take 
order 200 fb~^ to make a ha observation of a Mh — 150 GeV Higgs. 

5.6 Discussion 

The results summarized in Table 5.2 show that it is possible to isolate a virtually back- 
ground free qq qqH, H ^ tt signal at the LHC, with sufficiently large counting rate 
to obtain a 4 — 5(7 signal with a modest 60 fb~^ of data at low luminosity over most of 
the mass range. To reach 5(7 for Mh — 150 GeV would require about 200 fb~^. The 
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expected purity of the signal is demonstrated in Fig. 5.4, where the reconstructed rr in- 
variant mass distribution for a SM Higgs boson of mass 120 GcV is shown, together with 
the various backgrounds, after apphcation of aU cuts discussed in the previous Sections. 
This purity is made possible because the weak boson fusion process, together with the 
H T'^T~ f^hadrons^^T decay, provides a complex signal, with a multitude of 
characteristics which distinguish it from the various backgrounds. 

Additional cuts beyond forward tagging as discussed in Section 2.2 are specific to the 
if — > TT channel, with one r decaying leptonically and the other one decaying hadroni- 
cally. Crucial are charged lepton isolation and efficient identification of the hadronically 
decaying r, which are needed for the suppression of heavy quark backgrounds and non-r 
hadronic jets. This part of the analysis I have adapted from Ref . [59] , which, however, 
was performed for A, H ^ tt events from gluon fusion, i.e., without requiring two 
additional forward tagging jets. A more detailed assessment of lepton isolation and 
hadronic r identification in the present context is beyond the scope of the present work 
and should be performed with a full detector simulation. 

The elimination of the Wj + jj reducible background depends highly upon the 
Jacobian peak in the transverse mass distribution of the W decay products. The other 
backgrounds and the Higgs signal typically produce rather small values of mT{i,!^T), 
below 30 GeV, and thus well below the peak in mT{W). 

Another distinguishing feature of real r decays are the reconstructed momentum 
fractions x^-, and x^-^ of the charged decay lepton and of the decay hadrons. Misidentified 
"r's" tend to produce unphysically large values for these momentum fractions and can 
thereby be eliminated to a substantial degree (see Fig. 5.3). The reconstruction of these 
T momentum fractions is possible since the r+r" pairs are typically being produced 
with sizable transverse momenta (see Fig. 5.6c). As a result back-to-back r+r" decay 
products are rare (see Fig. 5.6b) and this in turns allows the mass reconstruction of the 
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T-pair, which is crucial for the suppression of the main physics background, Z — > rr. 

I have not made full use of the differences between the Higgs signal and the various 
backgrounds in some of these distributions. Additional examples are shown in Figures 
5.5 and 5.6. Fig. 5.5 shows the pTi and 77^ distributions for the observable charged 
lepton, which will form an important part of the event trigger. As a result of the lepton 
isolation cut, the falloff is considerably steeper for the bbjj background than for 
the signal and the other backgrounds. Not much leeway is present in applying more 
stringent cuts, however, without losing a substantial fraction of the signal. One can 
also take advantage of the if]i distribution for the QCD Zjj background, which, at the 
final level of cuts, remains important in particular for small values of the Higgs boson 
mass. 

In addition to the lepton py, one may use the missing transverse momentum of 
the event, ^t, Fig- 5.6(a), which is exceptionally small for the hbjj background. In 
combination with a more stringent cut on the r pair opening angle, cos(0,-^), shown 
in Fig. 5.6(b) (where an even more striking distinction between the physics and the 
reducible processes is found), both the Wj + jj and bbjj backgrounds can be reduced 
even below the level discussed in Section 5.4. Such a strategy, however, may not increase 
the statistical significance of the signal. In fact I find that shghtly looser cuts, for 
example on the dijet invariant mass, rrijj, can somewhat increase the significance of the 
signal while reducing the signal-to-background ratio. These points demonstrate that I 
have not yet optimized the search strategy for H ^ tt decays. This might be possible 
by combining the information from all the distributions mentioned above in a neural-net 
analysis. It is premature at this stage, however, to perform such an analysis since the 
issues of T-identification or of suppression of heavy quark decays in a realistic detector 
need to be addressed simultaneously, for the specific processes considered here. 
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Beyond the utility of confirmation of The Higgs' existence, independent measure- 
ment of a Higgs-fermion coupUng will be another important reason to strive for observa- 
tion oi H ^ TT decays at the LHC. For such a measurement, via the analysis outlined 
in this paper, r-identification efficiencies, minijet veto probabilities etc. must be pre- 
cisely known. For calibration purposes, the presence of the Z ^ tt peak in Fig. 5.4 
will be of enormous benefit. The production rates of the QCD and EW Zjj events can 
be reliably predicted and, thus, the observation of the Z — > rr peak allows for a direct 
experimental assessment of the needed efficiencies, in a kinematic configuration which 
is very similar to the Higgs signal. 
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M^V- (GeV) 

Figure 5.4: Reconstructed r pair invariant mass distribution for the signal and back- 
grounds after the cuts of Eqs. (2.1-2.3, 5.6-5.12) and multiphcation of the Monte 
Carlo results by the overall detector efficiencies (tagging jet identification efficiency 
^tag — (0.86)^ — 0.74, times r ID or jet/6 rejection efficiencies (see Eqs. (5.1, 
5.2, 5.5)) and expected survival probabilities. The solid fine represents the sum of 
the signal and all backgrounds. Individual components are shown as histograms: the 
Hjj signal (solid), the irreducible QCD Zjj background (dashed), the irreducible EW 
Zjj background (dotted) , and the combined Wj -\- jj and bbjj reducible backgrounds 
(dash-dotted). 
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Figure 5.5: (a) Transverse momentum and (b) pseudorapidity distributions of the 
charged "r" decay lepton after the cuts of Eqs. (2.1-2.3, 5.6-5.12), for the uih — 120 GeV 
signal (sohd hne), and backgrounds: QCD Zjj production (dashed hne), EW Zjj events 
(dotted hne), Wj+jj events (dot-dashed hne), and bbjj production (dash-double dotted 
hne) . 
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Figure 5.6: Shape comparison of various distributions for the Higgs signal (solid line) 
and the backgrounds: QCD Zjj production (dashed hne), Wj+jj events (dot-dashed 
line), and bbjj production (dash-double dotted line). Shown are the (a) -^t, (b) cos(0t-t) 
and (c) transverse momentum distribution of the reconstructed rr system, after the cuts 
of Eqs. (2.1-2.3, 5.6-5.12). 
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Chapter 6 
Conclusions 

The wonderfully developed and accurate theory of the electromagnetic, weak and strong 
forces of nature that is the Standard Model is nearly fully verified. It remains to be 
determined that the posited mechanism of spontaneous electroweak symmetry breaking 
and fermion mass generation, the Higgs mechanism, is correct. It is possible that some 
dynamical mechanism is instead responsible, and our precision electroweak data cannot 
determine this. It is also possible that the Higgs mechanism is partially correct: perhaps 
it provides for the symmetry breaking but not the fermion masses, or the other way 
around. Regardless, and barring any near-future discovery of a Higgs scalar at either 
LEP or the Tevatron, the LHC will have the capability to find a Higgs regardless of its 
possible mass. 

Of course, we are most interested in the lower end of the mass range, since precision 
electroweak data suggest the Higgs lies there. I have presented in this dissertation three 
additional modes in which to search for a Standard Model (or Standard Model-like) 
Higgs, via weak boson fusion in association with two hard tagging jets, providing for an 
experimental signature different from that of the inclusive (gluon fusion) searches. These 
modes completely cover the mass range from 100-200 GeV, safely overlapping both the 
LEP and Tevatron searches at the lower end and the "gold-plated" H ZZ — ^ Ai 
mode at the upper end. Furthermore, and importantly, all modes I consider here allow 
for Higgs mass reconstruction. 
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The first decay mode, H — > 77, covering the mass range of about 110-145 GeV with 
40-50 fb~^ at low luminosity (10^^ cm^^ s^^) for a 5cr statistical significance (the mass 
range can be expanded with additional luminosity, see Fig. 6.1), nicely complements 
the same search in gluon fusion, which can be accomplished with about 20 fb~^ with 
the CMS detector. At higher machine luminosities there will be an additional efficiency 
factor due to event rejection from minimum bias events. As the machine luminosity 
profile over time is not known, I do not attempt to include any factor in my estimates. 
In the event that nature doesn't generate fermion masses via the Higgs mechanism, this 
mode in weak boson fusion becomes extremely important to ensure that the Higgs does 
not go undetected. Mass resolution for this mode will be on the order of a GeV. 

The second decay mode, H — > — > e'^/jr^'^T, is much more complicated a 
search channel but also highly unique in its event characteristics. Observation of this 
mode is possible for very low integrated luminosities, on the order of 2-10 fb~^, if the 
Higgs boson lies in the mass range between about 130 and 200 GeV. 5a observation can 
be pushed down to about 115 GeV with about 200 fb~^ at low luminosity, see Fig. 6.1. 
There is even room for improvement from my presentation, shouki a multivariate anal- 
ysis with full detector simulation be performed on the observables I discussed. 

The third decay mode, H — > t+t~ — > i^h'^'^T, can cover the Higgs mass range 110- 
140 GeV with reasonable integrated luminosity, about 60 fb"^ (~ 30 fb~^ for Mh ~ 
115 — 130 GeV), and the mass range can be extended to 150 GeV with order 200 fb~^ 
of data, see Fig. 6.1. This search is also challenging, but I have shown that it is feasible 
and desirable, as it would provide the first direct Higgs-fermion coupling measurement. 
It will also be an extremely important mode for investigating the possibility of a MSSM 
Higgs sector, as it (combined with the WBF H ^ jj mode) can cover the entire 
tan P — Ma parameter plane. If either of /i or if is not observed, the MSSM is in serious 
jeopardy. 
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In all three cases, I advocate taking advantage of an additional fundamental charac- 
teristic of QCD and EW processes. Color-singlct exchange in the t-channel, as encoun- 
tered in Higgs boson production by weak boson fusion (and in EW Vjj backgrounds), 
leads to additional soft jet activity which differs strikingly from that expected for the 
QCD backgrounds in both geometry and hardness: gluon radiation in QCD processes is 
typically both more central and harder than in WBF processes. I exploit this radiation, 
via a veto on events with central minijets of pr > 20 GeV, and estimate a typical 70% 
reduction in QCD backgrounds, and a 20 — 25% suppression in EW backgrounds, but 
only about a 10% — 15% loss of the WBF Higgs signal. 

Beyond the possibility of discovering the Higgs boson in the H WW mode, or 
confirmation of its existence in the others, measuring the cross sections in both weak 
boson and gluon fusion will be important both as a test of the Standard Model and as a 
search for new physics. For such measurements, via the analyses described here, minijet 
veto probabilities must be precisely known. For calibration purposes, one can analyze 
Zjj events at the LHC. The production rates of the QCD and EW Zjj events can be 
reliably predicted and, thus, the observation of the Z ^ ii peak allows for a direct 
experimental assessment of the minijet veto efficiencies, in a kinematic configuration 
very similar to the Higgs signal. 

Weak boson fusion at the LHC will be an exciting process to study, for a weakly 
coupled Higgs sector just as much as for strong interactions in the symmetry breaking 
sector of electroweak interactions. 
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Figure 6.1: Minimum integrated luminosity (fb~^) required to observe WBF Higgs pro- 
duction and subsequent decay in each of the three modes described in this dissertation: 



^ 77 (sohd); H W^*^W^*^ e±/i^j!^T (dashed); and H 



T ' T 



(dotted). Values quoted are for low machine luminosity (10^^ cm ^ s ^) and thus do 
not include additional efficiency factors due to minimum bias event rejection. 
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Appendix A 
Important Parameters 

I simulate pp collisions at the CERN LHC, \/s = 14 TeV. For all our numerical results 
I have chosen sin^^^^ = 0.2315, Mz = 91.19 GeV, and Gp = 1.16639 • 10"^ GeV'^, 
which translates into Mw — 79.94 GeV and a{Mz)~^ — 128.74 when using the tree- level 
relations between these input parameters. This value of Mw is somewhat lower than the 
current world average of ~ 80.35 GeV. However, this difference has negligible effects on 
all cross sections, e.g. the qq — > qqH signal cross section varies by about 0.5% between 
these two W mass values. The tree level relations between the input parameters are 
kept in order to guarantee electroweak gauge invariance of all amplitudes. For all QCD 
effects, the running of the strong coupling constant is evaluated at one-loop order, with 
as{Mz) — 0.118. We employ CTEQ4L parton distribution functions [65] throughout. 
Unless otherwise noted the factorization scale is chosen as = min(pr) of the defined 
jets, and the renormalization scale is chosen as the px of the final-state colored parton 
emitted from a QCD vertex. For processes with more than one QCD vertex, the overall 
as is taken to be the geometric mean of each vertex's individual asi^pr), thus taking 
into account both the relevant scale for the hard scattering process and the soft scale 
for additional gluon emission in the relevant 3-jet processes. Inclusive WBF Higgs 
production cross sections as a function of mass are shown in Table A.l, as well as the 
HDECAY [42] -corrected branching ratios to 77, WW, and rr that are used in this study. 
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Table A.l: Signal inclusive cross sections (pb) for Hjj events of various Higgs masses in 
pp collisions at ^/s — 14 TeV, and the HDECAY-corrected branching ratios to photon, 
W and tau pairs. 



Higgs mass (GeV) 




BR{H 77) 


BR{H - 


WW) 


BR{H tt) 


100 


4.8 


0.00154 






0.0792 


110 


4.4 


0.00190 


0.047 




0.0763 


120 


4.1 


0.00218 


0.139 




0.0681 


130 


3.8 


0.00221 


0.299 




0.0532 


140 


3.5 


0.00191 


0.499 




0.0347 


150 


3.3 


0.00134 


0.700 




0.0175 


160 


3.1 


0.00051 


0.924 




0.0033 


170 


2.9 




0.967 






180 


2.7 




0.936 






190 


2.5 




0.779 






200 


2.4 




0.737 
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Appendix B 
Detector Resolution 



The QCD processes discussed in the previous Chapters lead to steeply falling jet trans- 
verse momentum distributions. As a result, finite detector resolution can have a sizable 
effect on cross sections. These resolution effects are taken into account via Gaussian 
smearing of the energies of jets/6's and charged leptons. All components of p'^ are 
smeared by the same factor derived from the energy smearing, i.e., — — p^. We 
use 

AE 5.2 0.4 , 

for jets (with individual terms added in quadrature), based on ATLAS expectations [41]. 
For charged leptons and photons I use 

AE _ 



= Z70 



(B.2) 



E 

which is quite conservative compared to CMS expectations [41]. 

In addition, finite detector resolution leads to fake missing-transverse-momentum in 
events with hard jets. An ATLAS analysis [66] showed that these effects are well pa- 
rameterized by a Gaussian distribution of the components of the fake missing transverse 
momentum vector, with resolution 



a{^.,^y)^0A6-JJ2ETMd: (B.3) 
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for each component. The coefficient 0.46 is vahd only in the low-luminosity scenario, 
which is sufficient as the searches in Chapters 4 and 5 do not require more than 30 fb~^. 
In my calculations, these fake missing transverse momentum vectors are added linearly 
to the neutrino momenta. 
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Appendix C 



Minijet Approximations 



C.l Introduction 

First introduced in Section 2.3, the minijet veto is a powerful suppression tool to ex- 
tract small WBF process signatures from large QCD backgrounds. While the necessary 
information on angular distributions and hardness of additional radiation is available in 
the 3-jet and tt + jets processes discussed in Chapters 3-5, one must either regulate or 
reinterpret these divergent cross sections. Here I will represent cross sections for 2- jet 
processes by (72. These calculations are completely perturbative and well-behaved, and 
will be regarded as inclusive cross sections for the respective processes under consider- 
ation. Cross sections for 3-jet processes, which are in general the respective 2-jet hard 
scattering process with one additional gluon emission, and all crossing related diagrams, 
will be represented by 



These cross sections are divergent for low pr of the additional gluon emission. For 
WBF processes, — 02 typically for PT,soft ~ 10 GeV, whereas for QCD processes, 
this occurs typically around pT,soft ~ 40 GeV for the LHC processes considered in this 
dissertation. 




(C.l) 
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In this Appendix I explain the technical details of both the truncated shower ap- 
proximation (TSA) [39], used to regulate the 3-jet calculations, and the exponentiation 
model, which alternatively reinterprets the meaning of the divergent cross sections (73. 



C.2 The Truncated Shower Approximation 

When several soft gluons are emitted in a hard scattering event their transverse mo- 
menta tend to cancel, leading to a regularization of the small-pr singularity (where 
the pt is the recoil of the hard scattering system) which is present when considering 
single-parton emission only. In the TSA these effects are simulated by replacing the 
tree-level three-jet differential cross section, daj^, with 

(iagTS^ = daj"^ (1 - e-^Ts/plsA^ . (C.2) 

Here the parameter Ptsa is chosen to correctly reproduce the tree-level two-jet cross 
section, (72, within the cuts specified for the comparison, typically the core forward 
tagging cuts of Eqs. (2.1-2.3), 

Pt^ > 20 GeV , \r]j\< 5.0 , ARjj > 0.7 , 

|77x|<2.5, ARjx>0-7. (C.3) 

Vj,min + 0.7 < 77x1,2 < ■nj,max - 0.7 , Tjj^ ■ r)j^ < . (C.4) 

A77to5, = |77,-, -77,,| >4.4, (C.5) 

where X is any observable Higgs decay product; plus any additional cuts on the tagging 
jets that alter the hardness of the underlying event, such as increased Ptj or Mjj cuts. 
I.e., Ptsa is fixed by the matching condition 

(72 = / -7^ dpT3 ■ (C.6) 

Jo dpT3 
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Ptsa values are typically ^10 GeV for WBF processes such as the Higgs signal, but 
much higher for QCD processes, characteristically 40-60 GeV. 

Using daj^^ as a model for additional jet activity one can make an estimate of the 
probability of observing at least one additional jet with p^^. > p^*° in the central region. 
These events can then be vetoed, as their angular distribution is determined correctly. 
However, in the TSA only one soft parton is generated, with a finite probability to be 
produced outside the veto region of Eq. (2.7). Therefore the veto probabihty will never 
reach 1, no matter how low PT,veto ma-Y be. 

At small values of pT,vcto one would expect to underestimate the veto probability for 
QCD processes because the TSA does not take into account multiple parton emission. 
Because of this our TSA veto probability estimates may be regarded as conservative. 



C.3 The Exponentiation Model 

In the soft region gluon emission dominates, and one may assume that this soft-gluon 
radiation approximately exponentiates, i.e., the probability P„ for observing n soft jets 
in the veto region is given by a Poisson distribution, 

= ^ e-^ , (C.7) 

with 




where the unregularized three-parton cross section is integrated over the veto region 
of Eq. (2.6) and then normahzed to the two-jet cross section, (72. A rough estimate of 
multiple emission effects is thus provided by using 

PexpiPT,yeto) = 1 - Pq = 1 - e-"^^'— ) (C.9) 

for the veto probability. 



96 



Within the exponentiation model, n — ^ represents the average multiphcity of 
minijets in the central region, between the two tagging jets. Even if the exponentiation 
model is of only limited accuracy, the ratio of three- to two-jet tree-level cross sections 
gives the best perturbative estimate available of the minijet activity in hard scattering 
events. One finds that the average minijet multiplicity depends strongly on the hardness 
of the underlying event. This is why I make estimates of veto probability only at the 
most stringent level of cuts on the hard scattering process. 



C.4 Application 

To employ a minijet veto in this study, I must first establish an algorithm for selecting 
tagging jets in the three- jet simulations such that events are selected that reflect the 
phase space region of the hard scattering of the analogous two-jet events. Only when 
this is ensured will the matching condition of Eq. (C.6) or the multiplicity definition, 
Eq. (C.8), make sense. In our previous studies of the minijet veto technique [35,37], the 
matching condition (or calculation of n) were performed without enforcement of the 
forward tagging cuts of Eqs. (C.4,C.5), even though tagging jet candidates were chosen 
for the purpose of identifying the veto candidate; tagging jet candidates were selected 
as the two most energetic [37] or two highest- pT [35] defined jets {pt > 20 GeV), in 
opposite detector hemispheres. However, without the additional forward tagging cuts, 
this will overestimate the effectiveness of the veto, especially for QCD processes. Taking 
QCD Zjjj production as an example, with Z — > i'^i" , the tagging candidate selection 
will frequently allow low-invariant-mass hard scattering events with gluon radiation 
forward of the quark jets to contribute to cxs, where one of the tagging jets is the gluon 
instead of the hard quarks of the LO process. Because the pt spectrum of the central 
quark jets, one of which is then identified as a veto candidate, can be very hard, on the 
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order of 40 — 60 GeV, the minijet "gluon radiation" spectrum is shifted toward higher 

Pt values. 

A more reahstic estimate of the minijet pt spectrum is obtained by applying the 
matching condition (or calculating n) only in the phase space region where a comparison 
of signal and background will take place: after all acceptance cuts, determined at the 
two-jet level, have been imposed; especially the forward tagging cuts of Eqs. (C.4,C.5), 
which significantly alter the hardness of the event selection. 

Once the full level of cuts for a given search scenario are imposed, one may examine 
different tagging jet selection algorithms to optimize the veto. Ideally, the outgoing 
quarks would always be selected, so that the additional gluon radiation is always the veto 
candidate. In practice, this is impossible, but for the Higgs signal various algorithms can 
achieve "proper" quark tagging with about 75% efficiency, a high success rate. Briefly, 
these might be the two highest-pr jets, or the two jets closest to the reconstructed 
Higgs. Most algorithms have very little difference from each other in the case of the 
WBF signature. Thus, I choose an algorithm that allows more suppression of the QCD 
backgrounds. The final algorithm I chose is to select the highest-pr jet as the first 
tagging jet, since it will almost always be part of the hard scattering, and then choose 
the other tagging jet such that the event is more likely to pass the forward tagging cuts: 
look for jets with pt > 20 GeV in the opposite hemisphere, such that the candidate 
Higgs decay products are between the tagging jets, satisfying Eq. (C.4). This performs 
somewhat superior to merely choosing the two highest-p^ jets. 

Also in contrast to our previous studies [35,37], the veto candidates are defined jets 
{Pt > 20 GeV) anywhere between the tagging jets, 

Vtag,min ^ Vj,veto ^ Vtag,max ■ 

(C.IO) 

Previously, the veto candidate also had to be at least 0.7 units of rapidity away from the 
tagging jets, but the choice of Eq. (C.IO) allows for more suppression of the backgrounds 
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than the more restrictive selection. 

Once the tagging selection algorithm is established, 1 need an estimate of the veto 
survival probability, Psurv, for the Higgs signal. The WBF Higgs processes including 
additional gluon emission were first discussed in Section 2.3, and the determination of 
Ptsa and estimation of Psurv with the TSA method were calculated via pp Hjjj 
production with decays to both if — > 77 and H ^ rr with the full level of acceptance 
cuts for each decay mode as discussed in their respective Chapters. This includes the 
forward tagging acceptance cuts of Eqs. (C.3-C.5) in both cases. For if — > 77 decays, 
the additional cuts are those of Eq. (3.3): 

PT,,i,,) > 40,20 GeV, 

PT,(i,,) > 50, 25 GeV. (C.ll) 

In this case, I calculate n — 0.067 and Ptsa — 5.4 GeV. This translates to Pgurv — 0-94 
for the exponentiation model and Pgurv — 0.88 for the TSA. Once Ptsa is set, I may 
switch to a tagging algorithm that selects the two defined jets closest in rapidity to 
the reconstructed Higgs, and find instead Psurv = 0.86. The small differences there are 
reassuring and due to slight changes in the allowable phase space of events. Considering 
the nature of the approximations in the exponentiation model and the TSA, it is also 
highly reassuring that their respective results agree to such a degree. For the purposes of 
this study I take the average of the results, Psurv — 0.89. These results are independent 
of the Higgs mass within Monte Carlo errors. 

For H ^ TT decays, the additional cuts beyond forward tagging are those of 
Eqs. (5.6,5.8-5.12): 

Pt,,^ ,, > 40,20 GeV, (C.12) 
PT.,., > 20 GeV , pT.,,„, > 40 GeV , (C.13) 

COSdrr > -0.9, Xri^^^ > 0, (C.14) 
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mjj >lTeV, (C.15) 

mT(^,j^T) < 30 GeV, (C.16) 

Xri < 0.75, Xr^ <1, (C.17) 

Here I find n — 0.128 and Ptsa — 5.5 GeV, which gives Pgurv — 0.88 for the expo- 
nentiation model, Pgurv — 0.87 for the TSA with pytagging as described above, and 
Psurv = 0.84 for the TSA with "?7-tagging" as before. The average is Pgurv = 0.87. That 
the survival probability is slightly smaller here than for the — > 77 case is probably 
due to the increased hardness of the event via the cut Mjj^tags > 1 TeV: the average 
maximum tagging jet px is ~ 100 GeV for H — > 77 but ~ 140 GeV for H — > rr; 
and similarly for the minimum tagging jet pt- By pt balancing, the additional gluon 
radiation in the Mjj^tags > 1 TeV events will be somewhat harder on average, leading 
to a slightly increased chance of being vetoed. 

Because these results arc so similar, and the Mjj cut is softer for the case H WW, 
there I apply the average value of Pgurv = 0.89 as found for the case of decays to photons. 

I also need to estabhsh veto survival probabilities for the backgrounds. It is not 
necessary to do this for all backgrounds, only for characteristic classes, e.g. QCD and 
EW production of weak bosons in association with two tagging jets, ti + jets, etc. 

For QCD and EW production of weak bosons in association with jets I examined 
QCD and EW Zjj{j) production with subsequent decay Z — > rr, as discussed in 
Sections 4.2.5 & 4.2.6, with subsequent tau decay to e^fj,^ or i^h"^, reflecting the 
backgrounds to H ^ WW and H — > rr, respectively. For tau decays to e^/^^ all the 
cuts of Eqs. (4.9-4.14) were used in addition to forward tagging: 

PTi > 20 GeV, (C.18) 

mjj > 650 GeV , (C.19) 
0e^<lO5°, cos^eM>0.2, AR,^ < 2.2 , (C.20) 
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< 110 GeV , pTe,^ < 120 GeV , (C.21) 

veto if > and niz - 25 GeV < m^^ < + 25 GeV . (C.22) 

For the QCD processes I find n = 1.46 and Ptsa = 43.5 GeV, yielding the estimates 
Psurv,n = 0.23, Psurv,TSAipT) = 0.34 and Psurv,TSA{r,) = 0.30. The average is Pguw = 0.29. 
For the EW processes I find h — 0.25 and Ptsa — 11-2, translating to Psurv,n — 
0.78, Psurv,TSA{pT) = 0.76 and Psurv,TSA(ri) = 0.71. The average is Pgurv = 0.75. The 
intermediate value of Psurv for the EW Zjj case reflects the partial bremsstrahlung 
nature of the hard scattering process, which can allow radiation back into the central 
region. 

For the decays Z ^ tt ^ i^h"^, 1 use the forward tagging cuts and those of 
Eqs (C.12-C.17). For the QCD processes I find n = 1.43 and pxsA = 41.8, yielding 
the estimates Psurv,n = 0.24, Psurv,TSA{pT) = 0.32 and Psurv,TSA{ri) = 0.28. The average 
is Psurv — 0.28, essentially the same as that found for the cuts of the H WW 
scenario above. For the EW processes I find n = 0.16 and Ptsa = 6.7, translating to 
Psurv,n = 0.85, Psurv,TSA{pT) = 0.82 and Psurv,TSA{r)) = 0.78. The average is Pgurv = 0.82. 
This survival probability is slightly higher than that for Z ^ tt ^ e^/i^, which is due 
to the fact that while the average px distributions of the jets is the same for both 
cases, the higher Mjj^tags requirement has moved event selection to the region where the 
bremsstrahlung component of Z production is much less, further reducing the fraction 
of minijets expected in the central region. 

1 use the above values for all QCD & EW backgrounds consisting of one or two 
weak bosons plus two tagging jets, for the respective cuts considered, and also the 
QCD backgrounds Wj + jj and hb + jj. The last is not quite the same, but does 
share hardness characteristics; a separate determination is not possible since the matrix 
elements for bb + jjj do not yet exist. For the QCD jjjj backgrounds (including DPS) 
I use Psurv — 0.30, extrapolating the other results for a lower value of Mjj^tags] and for 
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the EW 77 J J background I use Psurv = 0.75. The latter value probably underestimates 
the effectiveness of the veto as it is a yet lower Mjj^tags case. Thus, 1 am conservative 
in estimating the total background rates in the analysis of Chapter 3. 

Finally, I must consider tt + jets processes separately, as the typical hardness of 
the underlying event is not necessarily similar to any of the above processes. For these 
processes it is simpler in practice to reinterpret the divergent higher-order cross sections 
in the context of the exponentiation model, rather than use the TSA, which would need 
to be redefined in light of the additional jets. The survival probabihty for tt events 
was determined using the ttj processes where the jet is allowed to be soft and the 
two b jets from top decay were identified as the tagging jets. Similarly, the survival 
probability for ttj events was determined using the ttjj processes, where one b jets is 
identified as a tagging jet and one light quark or gluon jet is identified as the other 
tagging jet; one light quark or gluon jet is then allowed to be soft. I found a veto 
survival probabihty of Pgurv — 46% for tt events and Pgurv — 12% for ttj events. Both 
of these results disagree with our other estimates of Pgurv for QCD processes. This may 
be understandable for tt events, as at tree level this component does no contain any 
t-channel gluon exchange processes, which all of the other QCD backgrounds do. I also 
observed that the additional radiation in tt events typically falls outside the central 
gap. I did not explore this any further as the tt component is negligible. That the value 
of Pgurv found for tij events with the exponentiation model is so much smaller than 
that for other QCD backgrounds within the TSA may be understood for two reasons: 
the exponentiation method always gives a lower value for Pgurv than the TSA for QCD 
processes; and Ref. [48] have shown that the off-shell top contribution to the tt + jets 
background is not negligible, but I do not include it here. As these two issues are not yet 
fully explored I prefer to remain conservative and apply the value Pgurv — 0.29 for the 
tt + jets backgrounds. Thus, it is possible that I am overestimating their contribution. 
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Table C.l: Summary of veto survival probabilities for = 20 GeV used in Chap- 
ters 3-5. 
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tt 


ttj, 


QCD 
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QCD 


QCD 


DPS 








tijj 


V{V)jj 


V{V)jj 


Wjjj 


bbjj 


1133 


1133 


0.89 






0.30 


0.75 






0.30 




0.89 


0.46 


0.29 


0.29 


0.75 








rrjj 


0.87 






0.28 


0.80 


0.28 


0.28 
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Appendix D 
Tau Decay and Reconstruction 



D.l Tau Decays 



This analysis critically employs transverse momentum cuts on the charged r-decay prod- 
ucts and, hence, some care must be taken to ensure realistic momentum distributions. 
Because of its small mass, I simulate the r decays in the coUinear approximation. The 
energy fraction z of the charged decay lepton in — > I^viVt is generated according to 
the decay distribution 

1 dVo 1 



1-^ 



(5 + 5^-4^^) + Xr(l + -z-8z') . (D.l) 



dz 3 

Here Xt denotes the chirality of the decaying r (which, for a negative helicity t~ or 
positive helicity r^, is given by = — 1 in the relativistic limit). Similarly the pion 
spectrum for 'k^v^ decays is given by 

Decay distributions for r — > pv^- and r — > axv^ are taken from Ref. [58]. I add the decay 
distributions from the various hadronic decay modes according to their branching ratios. 
The vector meson decays are simulated in the narrow width approximation, which is 
adequate for my purposes, where the energy fraction is that of the vector meson. The 
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decay of the Higgs scalar produces r's of opposite chirality, Xt+ — ~Xt- and this anti- 
correlation of the polarizations is taken into account in our study. For the Higgs 
signal, there is no correlation of this final state with Higgs production, as the scalar 
carries no spin information. For vector boson production and decay there is correlation, 
however, which I simulate exactly for the case of the QCD Z — > rr background and 
approximate in the EW Z — > rr simulation. 

D.2 Tau Reconstruction 

A collider experiment will measure the momenta of massless final-state particles in an 
event, as well as the missing transverse momentum. I may denote these quantities for 
two massless final-state particles as ki, k2, and ^t, respectively. For a hadron collider, 
the longitudinal component of missing momentum cannot be determined. For dual tau 
decays ki, k2 are the momenta of the charged particles in the decay, leptons or hadrons. 
I will ignore all other final-state particles in an event here and assume that the missing 

— * 

momentum comes entirely from the escaping neutrinos in the two tau decays. I 
also neglect the r mass, as the coUinearity condition already requires a r transverse 
momentum much larger than the mass. 

Once the conditions for a coUinear r decay approximation are satisfied, the only 
unknowns are the two fractions of parent r energy which each observable decay particle 
carries, which I denote by x^-. in the text and abbreviate as Xi,X2 here. If pi,p2 are 
the tau momenta before decay, then we may write conservation of momentum in the 
transverse plane: 

— * — * 

(PlT + P2t) ^ ^ + ^ ^ hr+hr+fr. (D.3) 
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The transverse momentum vectors are then related by 



(D.4) 



As long as the the decay products are not back-to-back, Eq. (D.4) gives two conditions 
for and provides the r momenta as Pi/xi. Splitting into x and y components, this 
may be rewrittien in matrix form. 



which may be inverted to the form 



/ 



- 1 



- 1 



(D.5) 



/ 




- 1 






V 


J_ 

X2 


- 1 



1 I k2jT^ - k2jTy ^ 

kiMy - kiM [ k.jTy - kijT. J 



(D.6) 



Once Xi, X2 have been solved, it is then simple to reconstruct the tau-pair invariant 
mass as follows: 



ml+^- = {pi+P2f = 2{pi-p2 + ml) 



ki ■ k2 2 

h m_ 

X1X2 



2ki ■ k2 

X1X2 



(D.7) 



This technique is discussed in more detail in Ref. [63]. 
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Appendix E 
Monte Carlo Programs 

E.l Program Structure 

A typical Monte Carlo calculation begins with initialization of the phase space inte- 
gration limits and other constants, such as fundamental parameters. A loop is then 
performed for a large number of iterations, until the sampling error becomes small, 
wherein a set of random numbers is generated which is translated into four-momenta 
for the initial- and final-state particles. If the phase space point is valid and the configu- 
ration passes the cuts imposed on the final-state configuration desired, then the matrix 
element is calculated and squared for each subprocess. Each subprocess is folded with 
the structure functions for the incoming hadrons before they are summed together to be 
multiplied by the phase space weight. At the end of the calculation, collected statistics 
may be output to histograms for various phase space variables. 

Our general Monte Carlo package structure is a suite of FORTRAN (77 or 90) pro- 
grams, where each program performs a specific task in the process described above. 
The main program which drives the integration is generally named after the process 
being calculated. For example, to calculate the Higgs signal, the suite consists of the 
following program files: 

qqhmain . f 
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koppln . f 
qqhqq.f 
tbv.f 
monaco . f 
ps . f 
smear . f 
cuts . f 
tautau. f 
f eiketaus . f 
calcqsqr .f 
in2sZH.f 
hist.f 
f unc . f 

All parameters that need to be easily changed are put in . inc files, which centralizes 
the location so a change has to be made only once. 

qqhmain begins by initializing all the necessary variables and routines, from cou- 
plings (koppln. f, qqhqq.f, tbv.f) to histograms (hist.f). A loop over adaptive 
Monte Carlo iterations (monaco. f (a modified version of VEGAS [67])) and phase space 
points (ps.f) involves smearing the momenta (smear. f) for the acceptance cuts al- 
gorithm (cuts.f), and any particle reconstruction algorithms may follow (tautau. f, 
faketaus.f). If the phase space point selected survives the acceptance cuts (some of 
which may be in the reconstruction procedures), then the true momenta are passed to 
the matrix elements routine (m2sZH.f), which first calculates (calcqsqr. f) the val- 
ues necessary for structure function calls (CERN pdfiib). The driver then submits the 
event information to the histogramming subroutine (hist.f), which compiles results 
with CERN HBOOK routines. All extra useful subroutines, such as lego plot calculations. 
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may go in an extra file (func.f). 

Historically wc kept the routines in separate files to keep recompiling times short, but 
the advent of much faster machines may mean that the ultimate optimization achieved 
is slightly less than optimal. All computer codes used in this research, both signal and 
background, are available from the author upon request. 

E.2 Matrix Element Generation 

Many present-day hadron collider signal searches require the calculation of extremely 
complicated background containing many suproccesses with hundreds of Feynman dia- 
grams; the number of diagrams grows roughly with A^!, the number of final-state parti- 
cles in the configuration. Many calculations would be intractable or simply take years 
to prepare without the aid of Feynman configuration generators such as GRACE [68], 
MADGRAPH [4], or COMPHEP [69]. I have found madgraph to be the most powerful 
program; it has made possible the calculation of backgrounds with up to 8 final-state 
particles in a reasonable amount of time. In using MADGRAPH, the user inputs a given 
subprocess and the program outputs FORTRAN code which calls the helicity amplitude 
calculation subroutines of the helas package to find an amplitude for each graph, and 
then performs [Alp for the entire subprocess, with appropriate averaging. The user 
must determine each subprocess that is relevant and generate the code for it separately. 
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Appendix F 
Implications for the MSSM 

F.l The Higgs Sector in the MSSM 

In the minimal supersymmetric extension of the SM (MSSM) the search strategy for a 
Higgs is less clear [70] than in the SM. The Higgs sector in the MSSM consists of two 
complex scalar doublets, one with hypercharge +1 and the other —1, which after ESB 
similar to that in Section 1.3 result in five physical Higgs states: two CP even mass 
eigenstates, h and H, a CP odd A, and charged Higgs bosons H^. In this Appendix I 
summarize the reach of WBF with subsequent decay to rr for observing the neutral CP 
even Higgses in the MSSM. I show that the WBF channels are most likely to produce 
significant h and/or H signals in the regions of MSSM parameter space left uncovered 
by the MSSM Higgs searches at LEP. 

Relevant features of the MSSM Higgs sector can be illustrated in a particularly sim- 
ple approximation [72]: including the leading contributions with respect io Gp and the 
top flavor Yukawa coupling, ht = mt/{vs(i). The qualitative features remain unchanged 
in a more detailed description. All our numerical evaluations make use of a renormal- 
ization group improved next-to- leading order calculation [42,73]. The inclusion of two 
loop effects is not expected to change the results dramatically [74] . Including the lead- 
ing contributions with respect to Gp and /it, the mass matrix for the neutral CP even 
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Higgs bosons is given by 



r2 



+ m| 



'/3 L 



/\/f2 

log^^^ + 



^SUSY 



-S/3C/3 

(1- 



















1/ 



12M|(y5y^ 



(F.l) 



Here denote sin /3, cos /3. The bottom Yukawa coupling as well as the higgsino mass 
parameter are neglected (/x <^ ^Iusy)- The orthogonal diagonalization of this mass 
matrix defines the CP even mixing angle a. Only three parameters govern the Higgs 
sector: the pseudo-scalar Higgs mass, m^, tan/3, and e, which describes the corrections 
arising from the supersymmetric top sector. For the scan of SUSY parameter space 
I will concentrate on two particular values of the trilinear mixing term, = and 
At — y/6Mc;ijsY^ which commonly are referred to as no mixing and maximal mixing. 

Varying the pseudoscalar Higgs boson mass, one finds saturation for very large and 
very small values of - either nih or niH approach a plateau: 



ml 



m 



H 



for niA — > 00 , 
for tua — > . 



(F.2) 



For large values of tan (3 these plateaus meet at mi ff m| + e. Smaller tan P val- 
ues decrease the asymptotic mass values and soften the transition region between the 
plateau behavior and the linear dependence of the scalar Higgs masses on m^. These 
effects are shown in Fig. F.l, where the variation of mh and mn with ttia is shown 
for tan/3 = 4, 30. The small tan/3 region will be constrained by the LEP2 analysis of 
Zh, ZH associated production, essentially imposing lower bounds on tan (3 if no signal 
is observed. ^ 



-'^ Although the search for MSSM Higgs bosons at the Tevatron is promising [75] I quote only the 
Zh, ZH analysis of LEP2 [76] which is complementary to the LHC processes under consideration. The 
LEP2 reach is estimated by scaling the current limits for C = 158 pb~^ and ^/s = 189 GeV [76] to 
C = 100 pb-i and ^/s = 200 GeV. 
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The theoretical upper hmit on the hght Higgs boson mass, to two loop order, depends 
predominantly on the mixing parameter At, the higgsino mass parameter jj, and the soft- 
breaking stop mass parameters, which I treat as being identical to a supersymmetry 
breaking mass scale: ttiq — mu — MgusY [73]. As shown in Fig. F.l, the plateau mass 
value hardly exceeds ~ 130 GeV, even for large values of tan/3, Msusy — 1 TeV, and 
maximal mixing [74]. Theoretical limits arising from the current LEP and Tevatron 
squark search as well as the expected results from Zh, ZH production at LEP2 assure 
that the lowest plateau masses are well separated from the Z mass peak. 

The production of the CP even Higgses in WBF is governed by the hWW, HWW 
couplings, which, compared to the SM case, are suppressed by factors sin(/? — a), cos(/?— 
a), respectively [77]. In the plateau region (large m^), the mixing angle approaches 
a — P — 7i/2, whereas in the ttih plateau region (small m^i) one finds a ~ — /3. This 
yields asymptotic MSSM coupling factors of unity for h production and | cos(2/3)| ^ 0.8 
for the H channel, assuming tan/3 ^3. As a result, the production cross section of 
the plateau states in WBF is essentially of SM strength. In Fig. F.l the SUSY cross 
sections for cx^qq qqh/H) are shown as a function of m^; these may be compared to 
an expected rate for a SM Higgs of about 0.35 fb. The WBF cross section is sizable 
mainly in the plateau regions, and here the h or H masses are in the interesting range 
where decays into 66 and t^t~ are expected to dominate. 

The h and H couplings to 66 and rr are also modified in the MSSM by trigonometric 
factors of {(3 — a), but the details are unimportant here. For effective production of h or 
H by WBF, we have sin^(/3 — a) 1 or cos^(/3 — a) 1, respectively. The coupling of 
the observable resonance to 66 and rr is essentially of SM strength in these cases. The 
SUSY factors for the top and charm couplings are suppressed at large tan [3 , however, 
which leads to 66 and rr branching ratios similar to or exceeding the SM values for a 
given mass. 
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Once additional off-diagonal contributions to the Higgs mass matrix are included, 
it is likely that the hrr and Htt couplings become highly suppressed for sin(2Q;) = 0, 
which can occur in a physical region of the MSSM parameter space, although where 
exactly this may occur is strongly dependent on approximations made in the perturba- 
tive expansion. If the observed Higgs sector turns out to be located in this parameter 
region, the vanishing coupling to 66, rr would render the total widths small. This can 
dramatically increase the /i/if — > 77 branching ratio, even though T{h/H 77) may 
be suppressed compared to the SM case. In this region, the total h/H ^ should be 
sufficient for the search of Chapter 3 to be applicable, thus covering any region where 
h/H ^ TT becomes impossible. Additional details may be found in Ref. [78]. 



F.2 Higgs Search in Weak Boson Fusion 

Using the SUSY factors of the last section for production cross sections and decay 
rates, one can directly translate the SM results into a discovery reach for SUSY Higgs 
bosons. The expected signal rates, aB{h/H tt) are shown in Fig. F.l. They can be 
compared to SM rates, within cuts, of aB[H tt) = 0.35 fb and aB{H 77) = 2 fb 
for niH — 120 GeV. Except for the small parameter region where the tt signal vanishes, 
and for very large values of (the decoupling limit), the 77 channel is not expected 
to be useful for the MSSM Higgs search in WBF. The tt signal, on the other hand, 
compares favorably with the SM expectation over wide regions of parameter space. The 
SUSY factors for the production process determine the structure oi a ■ B{h/H ^ tt). 
Apart from the typical flat behavior in the asymptotic plateau regions they strongly 
depend on f3, in particular in the transition region, where all three neutral Higgs bosons 
have similar masses and where mixing effects are most pronounced. 

Given the background rates determined in Ch. 5, which are of order 0.03 fb in a 
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20 GeV mass bin, except in the vicinity of the Z-peak, the expected significance of 
the h/H TT signal can be determined. 5a contours for an integrated luminosity 
of 70 fb~^ are shown in Figs. F.2,F.3, as a function of tan/? and m^. I include an 
additional efficiency factor of 0.8 for both the signal and all backgrounds due to pile-up 
for data taken beyond the first 30 fb~^. Here the significances are determined from the 
Poisson probabilities of background fluctuations. Weak boson fusion, followed by decay 
to T-pairs, provides for a highly significant signal of at least one of the CP even Higgses. 
Even in the low tan /3 region, where LEP2 would discover the light Higgs, the WBF 
process at the LHC will give additional information. Most interesting is the transition 
region, where both h and H may be light enough to be observed via their rr decay. 

I have shown that the production of CP even MSSM Higgs bosons in WBF and 
subsequent decay to r pairs gives a significant (> 5cr) signal at the LHC. This search, 
with ^ lOOfb ~^ of integrated luminosity, and supplemented by the search for /i/iJ — > 77 
in weak boson fusion, should cover the entire MSSM parameter space left after an 
unsuccessful LEP2 search, with a significant overlap of LEP2 and LHC search regions. 
The two CERN searches combined provide a no-lose strategy by themselves for seeing 
a MSSM Higgs. At the very least, the WBF measurements provide valuable additional 
information on Higgs couplings. 

The present analysis relies only on the typical mixing behavior of the CP even mass 
cigcnstates, and on the observability of a SM Higgs, of mass up to ~ 150 GeV, in WBF. 
This suggests that the search discussed here might also cover an extended Higgs sector 
as well as somewhat higher plateau masses, e.g. for very large squark soft-breaking 
mass parameters. Because decays into r pairs arc tied to the dominant decay channel 
of the intermediate mass range Higgs, h/H bb, the search for a rr signal in WBF is 
robust and expected to give a clear Higgs signal in a wide class of models. 
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Figure F.l: Variation of h/ H masses, couplings to W/ Z, MSSM/SM strength ratio and 
total signal rate, for the CP even MSSM Higgs bosons as a function of the pseudoscalar 
Higgs mass. The complementarity of the search for the hghter h (left column) and 
heavier H (right column) is shown for tan/? = 4, 30 (dashed, sohd lines). Other MSSM 
parameters are fixed to = 200 GeV, MgusY — 1 TeV, and maximal mixing. 
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Figure F.2: 5a discovery contours for h tt and H ^ tt in WBF at the LHC, with 
70 fb"^. An additional efficiency factor of 0.8 applied to the signal and all backgrounds 
due to pile-up is included beyond the ffist 30 fb^^. Also shown are the projected LEP2 
exclusion limits (see text). Results are shown for SUSY parameters as in Fig. F.l, for 
maximal trilinear mixing, At = ^/6MsusY■ 
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Figure F.3: Same as Fig. F.2, but for the case of no trilinear mixing, At = 0. 
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